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Preface 



The lattice parameter, or the density, stands on an equal 
footing with the temperature in the theoretical description of a 
solid, liquid, or gas. This close connection with the theory is the 
basis for the importance of experiments conducted under high 
pressure. As the ingenious and prolific work of P. W. Bridgman 
has demonstrated, a very wide class of experiments can in fact 
be performed under conditions of high pressure. In 1953, a new 
class of experiments was added when the author and E. M. 
Purcell showed that it was possible to perform nuclear resonance 
experiments to 10,000 bars. Later developments by W. M. 
Walsh and N. Bloernbergeii extended the frequency range into 
the microwave region. Very recently the author and J. D. Litster 
have shown that it is possible to carry out nuclear resonance 
experiments into the pressure region above 50,000 bars. Today 
the experimenter can perform microwave and radiofrequency 
magnetic resonance experiments which enable him to examine 
the microscopic details of electronic structure and atomic or 
molecular motion, while at the same time he varies both the 
fundamental parameters the pressure and the temperature 
over the range 1 bar < P < 10 5 bar, and 1.5K < T< 400K. 

This tract presents an account of experiments performed in 
this field in sufficient detail so that the reader can see both the 
power and the limitations of the method. In so doing we cover 
a wide range of physical phenomena : from the magnetic proper- 
ties of ferromagnetic metals to the self diffusion of molecules in 
liquids and gases. In each case the presentation has focused on 
the physical principles involved and the new information ob- 
tained by experiment. Generous literature references have been 
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provided to enable the reader to obtain further theoretical and 
experimental details from the original work. As an aid to the 
experimenter I have included in a short Appendix an account 
of some of the principal techniques that are used in magnetic 
resonance experiments at high pressure. In addition to an 
examination of the present state of the field, the author has 
endeavored both by organization of the material and by 
direct suggestion to point out still fruitful areas for 
experimentation. 

While the present material is devoted to magnetic resonance 
at high pressure, it is important to recognize that were we to ex- 
tend our attention beyond the radio frequency and microwave 
region we would find that optical spectroscopy of solids may be 
carried out to 100,000 bars, as the work of H. Drickamer has 
shown. Beyond the optical region, x-ray diffraction and 7-ray 
spectroscopy has been used successfully to study the properties 
of solids at high pressure. For example, we include in this tract 
a discussion of the effect of pressure on the energy of the 14.4-kev 
7-ray emitted by the Fe 57 nucleus in ferromagnetic iron. Today 
experiments on the properties of matter under high pressure 
have been conducted using electromagnetic radiation whose 
frequency ranges all the way from d. c. to 10 18 cps. 

In addition to this great breadth available in the electro- 
magnetic spectrum, the stress applied to the sample need not 
be confined solely to hydrostatic pressure. Important magnetic 
resonance experiments have been conducted by G. Feher, G. 
Watkins, R. G. Shulman, and others on solids subjected to uni- 
axial stress. However, in this short tract it will not be possible to 
include an account of their work. 

In the interests of clarity we include a note on the units most 
frequently employed in high pressure work the kilogram per 
cm 2 , the atmosphere, and the bar. These have the following rela- 
tions to one another and to other commonly used units: 1 kg/ 
cm 2 = 0.9678 atm = 0.9807 bar = 14.22 pounds/inch = 0.9807 
X 10 6 dynes/cm 2 . 
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Chapter I 

Solids at High Pressure 



1. Magnetic Fields in Solids 
1.1 Introduction 

Both nuclei with magnetic moments and unpaired electrons 
serve as sensitive probes of the magnetic fields in a solid, and 
both have been used to study these internal fields. 

In the nuclear case, a magnetic field // at the position of a 
nucleus interacts with the nuclear magnetic moment gN@ N I to 
produce a set of energy levels 



E(m) = -<72vftv//m 7 (1-1) 

where m/ is the magnetic quantum number corresponding to the 
nuclear spin I, and g N and N are respectively the nuclear g 
value, and the nuclear magneton (1) N = eh/2MC. By applying 
to the nuclear moments a perturbing radio-frequency field per- 
pendicular to //, resonant absorption of rf power may be ob- 
served when the frequency of the perturbing field corresponds 
to splittings for which Am/ = 1. Measurement of the resonance 
frequency v thus enables determination of the field H by the 
relation 

v = (g N (*N/h}H (1-2) 

In favorable cases the line width (A*>) of the nuclear absorption 
is very small compared to the resonance frequency, and hence a 
very accurate determination of H is possible. 

Apart from an externally applied field, the magnetic field at 
a nucleus is produced by both the orbital and the spin magnetic 

1 
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moments of the surrounding electrons. In solids the crystalline 
electric fields frequently quench the orbital angular momentum 
so that in first approximation we may confine our attention to 
the effects of the spin magnetic moments. Let T* be the position 
vector of the ith electron, r^ the position vector of the nucleus in 
question, and further let r iN - be the vector joining the two. Let 
HO be an externally applied magnetic field. The Haniiltoniun 
representing the interaction between the nuclear moment /^ v 
and the internal and external fields is given by 

/ ?5 \ 



where the sums are carried out over all the electrons and y* 
g@Si is the magnetic moment of the ilh electron (g = 2 and ft 
is the Bohr magneton). The second term in the cur ley brackets 
is the operator which represents the field set up at the nucleus by 
the magnetic moment of electrons situated at the nucleus. 
This is the hyperfine, or contact interaction, first obtained by 
Fermi (2,3). This term is usually much larger than the third 
term, the dipole-dipole term, which represents the magnetic/ 
field set up by electrons which are at points removed from the 
nuclear position r N . 

It is vital to remember that the internal fields produced by the 
electrons in a metal, paramagnetic, ferromagnet, or antiterro- 
magnet, fluctuate very rapidly compared to the Larmor fre- 
quency of the nuclei. Thus the nuclei respond only to the time- 
average of the field set up by the electrons, and the nuclear 
resonance frequency is shifted from g#ft N II Q /h by an amount 
^nt = g N N (Hini)/h. Thus 

v = QN&yrHt/h + v^t 

To determine (//mt) one must take a statistical average of the 
last two terms in equation (1-3) over the various stationary 
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states of the electronic system. The evaluation of this sta- 
tistical average depends in detail on the nature of the solid 
under study. Nevertheless, since the local fields are proportional 
to the magnetic moment ju of the electrons, we can expect (Hint) 
to be proportional to the magnetization M(T), which is a func- 
tion of the temperature (T). Hence 

(// int ) = A'M(T) (1-4) 

In paramagnets M is proportional to HQ. In feebly paramagnetic 
metals (Hi n t)/H is the Knight shift. In ferromagnets and anti- 
ferromagnets there is a spontaneous internal magnetization in 
zero external field, and the nuclear resonance can be observed 
in the absence of an external field. 

A' contains information on the unpaired spin probability 
density for electrons at the nucleus. (Two electrons which have 
the same spacial wave functions and differ only in spin states are 
called paired electrons. These have no net spin magnetic mo- 
ment, and hence make no contribution to 3C.) A ' also contains 
information on the dipole sum in equation (1-3). We therefore 
see that the shift 

*int = AM(T) (1-5) 



of the nuclear resonance frequency away from j> = 

presents us with information jointly on the electronic wave 

function and on the magnetization, M. 

It is clear from equation (1-5) that studies of the pressure- 
dependence of the internal fields, when combined with measure- 
ments of the compressibility (l/V)(dV/dP) T9 can provide in- 
formation on the volume-dependence of the coupling constant 
and the volume-dependence of the magnetization, since 



(d In ? int /d In V) = (d In iWdP)r/(V V)(5V/dP) r (l-6a) 
= (d In A/b In V) T + & In M/d In V) T (b) 

The volume-dependence of A contains information on the al- 
teration of the electronic wave functions as the atoms in the 



4 MAGNETIC RESONANCE AT HIGH PRESSURE 

solid are brought closer together. The volume-dependence of the 
magnetization, taken at different temperatures, contains in- 
formation on such quantities as the volume-dependence of the 
Curie, or Neel, temperature, and the volume-dependence of the 
saturation magnetization. 

A more subtle, but nonetheless important, application of 
measurements of the pressure-dependence of ^ int is in correcting 
measurements of the temperature-dependence of v\ n t for the 
effects of thermal expansion. In using equation (1-5) to obtain 
accurate information on the temperature-dependence of the 
magnetization it must be kept in mind that comparison with 
theory requires measurements of *>i nt to be made without change 
in the size of the unit cell. Since measurements of z>i n t are made 
at constant pressure, thermal expansion alters the size of the 
unit cell. This effect can be corrected using the pressure measure- 
ments. For simplicity, let us assume that we deal with a cubic 
crystal in which the dimensions of the unit cell and the positions 
of atoms inside it can be specified by a single parameter which, 
we shall take as the volume (V) of the unit cell. In this case the 
relation between constant-pressure and constant-volume meas- 
urements of the temperature-dependence of j>mt is given by 

(d In vint/bT) v = (d In v in t/t>T) P 

- (d In Pint/d In V) T (d In V/t>T) P (1-7) 

The last term in equation (1-7) is the effect of the thermal ex- 
pansion. The magnitude of this term depends on the size of the 
thermal expansion coefficient and the volume-dependence of 
p^t. This latter quantity can be obtained from the pressure- 
dependence of Vmt by equation (l-6a). Thus, if the compressi- 
bility (l/y)(dV/dP)r is known at various temperatures, 
measurements of the pressure-dependence of the internal field 
enables the determination of (d In Vmt/& In V) and hence the 
volume expansion correction in equation (1-7). Pressure meas- 
urements are necessary, as we see, in making the magnetic 
equivalent of the famous C P to C v corrections in the theory of 
specific heats. 
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The electrons, too, serve as probes for the internal fields in a 
solid. In Section 1.5 of this chapter we deal with the measure- 
ment of the pressure-dependence of the hyperfine interaction 
constant between the electrons and nucleus of paramagnetic 
ions substituted into various host lattices. 

1.2 Feebly Paramagnetic Metals (The Knight Shift) 

The magnetic field at the nucleus of an atom in a nonferro- 
magnetic metal can be determined from the Hamiltonian (equa- 
tion (1-3)), provided that the stationary states of the electrons 
are known. In a metal the electrons often can be unambiguously 
separated into two groups the core electrons, and the conduc- 
tion electrons. The core electrons nominally consist of electrons 
with spins paired off in spherically-symmetrical closed shells. 
Such electrons make no contribution to the field at the nucleus. 
The valence electrons are distributed throughout the metal in 
spin-up and spin-down Bloch states, and electrons occupy these 
states in accordance with the Fermi-Dirac statistics. 

It is necessary at the outset to estimate the correlation time, 
T C , for fluctuations in the field at the nucleus set up by the con- 
duction electrons. Because of the Fermi statistics it is necessary 
to consider only those electrons near the Fermi surface. The 
conduction electrons well below this energy have their spins 
completely paired. The electrons near the Ferrni surface have a 
velocity VF and jump from one unit cell to another in a time r c ~ 
a/v F ~ h/f where a is roughly the edge length of the unit cell, 
and f is the Fermi energy. It is this jumping from one atom to 
another which produces the fluctuation of the field at the nu- 
cleus. The correlation time, r c ~ 10~ 16 sec, is much shorter than 
the inverse of the Larmor frequency (l/v L ~ 10 ~ 7 sec) of the 
nucleus. As a result, we see at the outset that the field at the 
nucleus must be calculated by taking the statistical average of 
the term in curly brackets in the Hamiltonian (equation (1-3)). 

To carry out the statistical averaging, we introduce the 
density of conduction band states, n(e), around the energy c, and 
the Fermi f unction, / (e - f/kT). In the presence of a magnetic 
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field it is convenient to divide the density of states into moment- 
up and moment-down (n + (e) and n~(e)) parts. It is also neces- 
sary to include in the argument of the Fermi function the 
potential energy ##o of the electrons. If we now define N + 
and N~ as the total number of electrons with moments respec- 
tively parallel and antiparallel to the applied field, then the 
relation 



determines the Fermi level f . By expanding the Fermi function 
in powers of (3Ho/kT it is seen that to first order the Fermi level 
is unaffected by the application of a magnetic field. To carry out 
the statistical average of the magnetic field at the nucleus (H N ) 
we observe that for cubic metals the dipole-dipole terms in equa- 
tion (1-3) give no contribution. The hyperfine interaction term 
gives a contribution only for the conduction electrons, and we 
obtain 

- H 

N g 



/ + 8//o - r\i . 

- rc~(e)/o ( -- - 1J de 



n m 

(1-9) 



The quantity in square brackets is zero, except for a narrow 
range of the order of kT near the Fermi energy f. Thus if 
\&(r N , c)| 2 is a slowly- vary ing function of the energy e, we may 
remove it from the integral and set c = f in its argument. The 
terms that remain inside the integral give exactly the magnetic 
moment M 
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M 



- ..-(<)/. (' d, (1-10) 

Thus we may write 

II N = Ho + ^- x*|*(r w , *)| 2 #o = #o(l + K) (1-11) 

t> 

The Knight shift K may be written as 

K = (Bir/SKxpOPiO (1-12) 



where all the quantities entering into equation (1-12) have been 
normalized per atom. %P is the magnetic susceptibility per unit 
volume, ft is the volume per atom, and PF = \&(rn, f ) | 2 where ^ 
is normalized over the unit cell. The quantity % P 12 is the para- 
magnetic susceptibility per atom, given in the free-electron 
model by the Pauli expression 

(1-13) 



where p = number of valence electrons per atom. It will be ob- 
served that the contribution to the resonance frequency which 
comes from the internal field is exactly of the form of equation 
(1-5), viz. 



where P F plays the role of the coupling constant and XP^#O is 
the magnetization, M, per atom. 

Before turning to the measurement of the pressure-dependence 
of K, it is well that we consider an important effect that occurs 
again and again in the study of the hyperfine interaction, namely 
core polarization (4,5). We have heretofore explicitly assumed 
that the core electrons are distributed in orbitals in which the 
spins are exactly paired off, i.e. for each electron in a spin-up 
state there is another in an identical spin-down state. This as- 
sumption will be correct if there are equal numbers of spin-up 
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and spin-down electrons in the system under consideration. 
However, in magnetic materials this is not the case. As Slater 
has emphasized (6) the existence of a net spin in the system 
means that a spin-up electron (in our case a spin-up core elec- 
tron) will see a different exchange potential than will a spin- 
down electron. The result is that spin-up and spin-down core 
electrons will have slightly different orbitals and, in particular, 
the spin-up and spin-down probability densities at the nucleus 
will not be identical. Thus the core electrons are polarized by the 
exchange interaction with the magnetic electrons and con- 
tribute to the total unpaired spin density at the nucleus. No 
explicit (7) calculations have yet been made for the effects of 
core polarization on the Knight shift in metals. Nonetheless, 
we must keep in mind that PF as deduced from the Knight shift 
does contain contributions from core electrons polarized by ex- 
change interaction with magnetic electrons at the Fermi level f . 
The pressure-dependence of K has been measured (8) in 
lithium, sodium, rubidium, cesium, copper, and aluminum at 
room temperature over the pressure range 1 to 10,000 kg/cm 2 . 
The measurements consist of the determination of the change 
in nuclear resonance frequency v when hydrostatic pressure is 
applied to a sample placed in a constant external field, HQ. If 
v(P) and *>(0) are respectively the nuclear resonance frequencies 
at pressure P and at 1 atmosphere it follows from equations 
(1-2) and (1-11) that 

K(P) - K(0) = ^"^ (1 + K(0)) (1-15) 

Thus, measurement of the pressure-dependence of K requires the 
determination of the change in nuclear resonance frequency on 
applying the pressure, and the values of K at atmospheric pres- 
sure. The latter quantities are well known (9). The changes in v 
however are very small in the cases of lithium and sodium. At 
10,000 kg/cm 2 , the change in v is only ~30 cps for lithium and 
110 cps for sodium. Hence, stringent requirements are placed 
on stabilization of the external field and the frequency stability 
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Figure 1. The change in the Na 23 nuclear resonance frequency as a 
function of pressure at 24.4C. The accuracy of the determination of the 
Na frequency is 1 cps at each pressure. 

of the spectrometer. To meet these requirements, a system, then 
unprecedented in accuracy, was developed which enabled 
stabilization of both spectrometer frequency and magnetic field 
to a few parts in 10 7 over long periods of time (8 to 10 hours). 
The results for sodium are shown in Figure 1. It will be noted 
that the center of the resonance line is determined to within 1 
cps at each pressure. When the results for the pressure de- 
pendence of K are combined with Bridgman's measurements of 
the volume compression of the solid we obtain the volume- 
dependence of K . 

The primary object of the measurement of the volume-de- 
pendence of K was the determination of the volume-dependence 
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of Pp. This may be obtained from the volume-dependence of the 
Knight shift, provided that the volume-dependence of XP is 
known. Clearly it would be most desirable to make an inde- 
pendent experimental measurement of the volume-dependence 
of the Pauli susceptibility. Unfortunately, however, the experi- 
mental accuracy now available for measurements (10) of XP 
are too low to permit determination of its volume-dependence. 
Instead, the collective electron picture of Pines (11), which is 
successful in predicting quantitatively XP at atmospheric pres- 
sure, was used to calculate the volume-dependence of XP^- 
Combining this with the measurements of K vs. volume, the 
volume-dependence of P F was deduced and compared with the 
theoretical predictions of W. Kohn and F. Ham and H. Brooks. 
The quantum defect method developed by Brooks and Ham 
(12) permits calculation of the volume-dependence of PF for all 
the alkali metals without explicit construction of the potential 
in which the conduction electrons move. This theory gives very 
good results (8) for the tolume-dependence of P F for sodium and 
rubidium. Cesium shows a much stronger volume-dependence 
than is predicted by the theory. Because of convergence difficul- 
ties the theory gives no prediction for lithium. These results are 
summarized in Table 1-1. 

The temperature-dependence of K is also of interest. McGar- 
vey and Gutowsky (13) measured K vs. the temperature (7 1 ) for 
Li, Na, Rb, Cs. They believed that the temperature-dependence 
was due entirely to the thermal expansion: that K was not an ex- 
Table 1-1 

The Fractional Change in the Knight Shift and P F on Changing the 
Volume of the Solid by the Fraction 



Metal AV/V AK/K AP F /P P (theoretical) 



Lithium 


-0.08 


-0 0116 


-0.005 




Sodium 


-0.13 


-0.0131 


+0.072 


+0.062 


Rubidium 


-0.23 


+0.068 


+0.22 


+0.21 


Cesium 


-0.24 


+ 0.415 


+0.57 


+0.34 
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plicit function of the temperature. In terms of equation (1-7), in 
which Pi n t is replaced by K, they believed that (d In K/dT) v = 0. 
The availability of the pressure measurements discussed above 
enables the determination of (d In K/d In V)T> When this is 
combined with the known thermal expansion (d In V/dT) P 
and compared to (5 In K/bT) P , it is found that K is an explicit 
function of the temperature at constant volume. It was pro- 
posed (8,8a) that this effect is due to the modulation of P F by 
the lattice vibrations. These vibrations cause a fluctuation 
of the volume of the unit cell, and a concomitant fluctuation 
in the value of PF. The susceptibility per atom %p^ does not 
follow the very rapid fluctuation because the time required for a 
spin flip is long compared to the characteristic period of the 
lattice vibrations. The nucleus therefore responds to the time- 
average value of P F , which is 



(1-16) 



The magnitude and temperature-dependence of (12 ) 2 can be 
computed using the Debye model of lattice vibrations, while 
d 2 PF/d2 2 can be obtained from the volume-dependence of PF as 
deduced from the pressure-dependence of K. Calculations (8) 
using this model are in good agreement with the experimental 
results. 

1.3 Antiferromagnets 

Manganous fluoride, MnF2, is a nearly ionic crystal. On form- 
ing the crystal lattice, in first approximation, the manganese 
atoms become doubly-ionized, magnetic, Mn 2 + ions with the 
configuration 6 S 5 /2, while the fluorine atoms become singly- 
charged nonmagnetic F~ ions in the closed shell configuration 
1 S Q . The magnetic moments on the Mn 2+ sites are responsible for 
the paramagnetism of MnF 2 at temperatures above the Neel 
temperature, T N . Below T N the super exchange interaction and 
anisotropy fields produce an antiferromagnetic ordering of the 
Mn 2+ spins. In Figure 2 we see the arrangement of nearest- 
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Figure 2. The arrangement of nearest-neighbor Mii 2 + ions around a 
fluorine nucleus in MnF2. The arrows represent the magnitude and relative 
orientation of the time-averaged magnetic moment of the Mn 2 + ions. A 
second-type fluorine site sees the same configuration as above, except that 
the direction of the Mn 2 + spins are reversed. 

neighbor Mn 2 + spins around each fluorine site. The nucleus 
under study is the F 19 nucleus. In this crystal there are two types 
of fluorine sites. In the antiferromagnetic state, one type of 
fluorine sees the spin arrangement shown in Figure 2, while the 
other type sees a reversed orientation of manganese spins. 

To examine the magnetic effect of the near-neighbor Mn 2+ 
spins on the F 19 nucleus we employ the Heitler-London, or 
atomic orbital, description of the wave functions in the crystal. 
If we neglect entirely the effect of overlap between the fluorine 
and manganese orbitals, then we may expect that the field (H) 
that the fluorine nucleus sees is due entirely to the magnetic 
dipole field of the Mn 2 + ions. At T = K, in the case of com- 
plete antiparallel alignment, this would produce a magnetic 
field of 12,500 gauss, 12,800 of which is due to the nearest-neigh- 
bor ions and 300 of which is contributed by the further neigh- 
bors. V. Jaccarino and R. G. Shulman (14), who were the first 
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to observe the F 19 nuclear resonance in MnF 2 , found that at 
T ~ K, the magnetic field at the fluorine site was actually 
40,470 gauss. They proposed (15) that this large field was due 
to an indirect hyperfine interaction between the spin of the Mn 2+ 
ions and the F 19 nucleus, and suggested the possibility that the 
observed hyperfine interaction was the result of covalent bond- 
ing between the fluorine ion and its Mn 2+ neighbors. The 
covalent bonding requires a transfer of fluorine electrons onto 
the manganese neighbors. The transferred electrons must have 
their spins antiparallel to those on the Mn 2+ neighbors, since 
each of these already has its full complement of 5 parallel 
spins. The unpaired spins remaining on the F~ ion then interact 
with the F 19 nucleus through the hyperfine interaction. Despite 
this suggestion the origin of the Mn-F hyperfine interaction 
was not clearly established at this stage in the development of 
the subject. 

Turning now to the situation at temperatures above K, we 
must examine the correlation time (r c ) for fluctuations of the 
Mn 2+ spins. This correlation time is determined by the Mn-Mn 
superexchange interaction which couples each pair of Mn 2 + 
ions in accordance with 3Ci 2 = JSr 82. Fluctuations between the 
two degenerate states (1 T), (2 i) and (1 |), (2 t) take place 
with a frequency w e ~ J/h. Estimating the exchange integral J 
from the Neel temperature kT N ~ J, we have o? 6 ~ 7V/h. 
Using T N ~ 68 K (16) we find co, ~ 8 X 10 12 /sec and r c = 
27r/u e ~ .8 X 10~ 12 sec. Clearly this is much shorter than the 
inverse of the Larmor frequency of the fluorine nucleus, so 
again we see that the nucleus responds only to the time-average 
spin orientation of the neighbor manganese spins. This time- 
average spin orientation, in turn, is proportional to the sublat- 
tice magnetization, M(T). Thus we may write that the resonance 
frequency, i>, is related to the sublattice magnetization by 

v = AM(T) (1-17) 

which of course is identical in form with equation (1-5). The 
constant of proportionality, A, in equation (1-17) includes both 
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the hyperfine coupling constant and the dipole-dipole inter- 
action. 

At temperatures well below the Neel temperature, M(T) be- 
comes a constant independent of the temperature and very 
nearly independent of the volume. Thus, by measurements of 
the pressure-dependence of the resonance frequency v at very 
low temperatures and in zero external field, one can determine 
the pressure-dependence of A, and hence the pressure-de- 
pendence of the hyperfine coupling constant. Furthermore, if 
the change in lattice dimensions produced by the pressure is 
known, one may determine how the hyperfine coupling constant 
depends on the distance between the manganese and fluorine 
ions. This experiment has been performed (17) at 4.2 K using 
liquid helium as a pressure-transmitting medium. The pressure 
range was limited to only 100 kg/cm 2 , but the accuracy with 
which the line center could be determined was sufficient to 
specify that 

(1/AX&A/&P) = (1.9 0.1) X 10~V(kg/cm 2 ) 

At higher temperatures, where M(T) is not a constant, the 
pressure-dependence of v arises not only out of the pressure- 
dependence of A, but also out of the pressure-dependence of the 
magnetization, the sublattice magnetization being a function of 
pressure through its dependence on the Neel temperature T N . 
We may therefore write 



(1-18) 

If we take M = M(T/T N ) as the functional form of the de- 
pendence of M on T and T N , it follows that 



(1-19) 
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Since (d In M/dT) TN ^ (d In ?/dT)p, we may obtain (dT N /dP) 
by measuring the pressure-dependence of v at a temperature at 
which the second term on the right-hand side of equation (1-19) 
is appreciable compared to the first. This measurement has been 
carried out at 35.7 K, yielding a determination of 

(dT N /dP) = (0.298 0.02) X 10~ 3 K/ (kg/cm 2 ) 

using TV ~ 68.5 K. As a check of self-consistency, we can 
measure (\./v)(bv/'bP) T at some other temperature and compare 
the result with the prediction of equation (1-19) obtained using 
the values of d In A/dP and d In T N /dP given above. At 20K 
equation (1-19) predicts (I/v)(b/bP) T = (2.2 =fc 0.1) X 10~V 
(kg/cm 2 ) whereas we found (2.21 0.04) X 10 ~ 6 / (kg/cm 2 ). 
This satisfactory agreement lends support to the use of equation 
(1-19). 

Recently (18) the pressure-dependence of the Neel tempera- 
ture in MnF 2 has been determined very accurately by making 
use of the fact that at the Neel temperature the F 19 nuclear 
resonance line width broadens so drastically that the line in- 
tensity falls beneath the noise. Any change in the temperature 
at which this disappearance takes place could be detected to 
within 2 millidegrees. By measuring the pressure-dependence 
of the disappearance temperature we have found 

(dT N /dP) = (0.302 .002) X 10~ 3 K/ (kg/cm 2 ) 

which is in excellent agreement with the results obtained above. 
With the aid of measurements of the compressibility of MnF 2 , 
these determinations of dT N /dP show how the superexchange 
interaction between Mn 2+ ions depends on the distances and 
angles which characterize the Mn-F-Mn geometry. No theoreti- 
cal calculations have yet been made to predict the dependence 
of the superexchange integral on the Mn-F-Mn distances and 
angles. 

On the other hand, a theory has been proposed (17,19,20) to 
account for the magnitude and pressure-dependence of the 
hyperfine interaction between manganese ions and the fluorine 
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nucleus. In essence, this theory proposes that the unpairing of 
the core electrons in 5 states around the fluorine nucleus arises 
out of Pauli, or exchange, repulsion between fluorine and man- 
ganese spins which have the same orientation. Because of 
orthogonality, this exchange repulsion docs not act on those 
fluorine electrons whose spin is antiparallel to the direction of 
the neighbor Mn 2+ spin. As a result, spin-up and spin-down s 
and p orbitals around the fluorine nucleus are slightly different. 
The situation described here is physically quite similar to the 
core polarization effect mentioned in connection with the 
Knight shift. In the present case, the core is polarized by elec- 
trons whose orbitals are centered on sites different from the 
nucleus experiencing the hyperfiiie interaction. In the more 
familiar example (4,5,7) the core is polarized by unpaired elec- 
trons centered on the nucleus in question. Furthermore, in the 
present case, it is possible (17) to calculate the hyperfine interac- 
tion produced by the Pauli repulsion kinematically without 
constructing the exchange potentials and solving the unre- 
stricted Hartree-Fock equations. 

Using the model of Pauli or exchange repulsion first proposed 
by Mukherji and Das (20), along with W. Marshall's estimate 
of the alteration of the Hartree-Fock free-ion wave functions 
appropriate to the ions in the solid, calculations of the magnitude 
and pressure-dependence of A were made which were in excel- 
lent agreement with observation (17). This agreement lent sup- 
port to the view that covalency effects were not important in 
impairing the s orbitals on the fluorine ion. However, subse- 
quent refinements of the theory (19,21) showed that a cross- 
term between the 25 and Is orbitals had been omitted in the 
original calculation (17) . This term reduces the magnitude of the 
theoretical prediction of the hyperfine coupling constant by a 
factor of two. This difficulty may be due to inadequate knowl- 
edge of the details of the 3d Mn 2+ orbitals, especially in the 
region near the fluorine ion. On the other hand, we must now 
realize that covalency may play a significant role in determining 
the hyperfine interaction. 
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1.4 Ferromagnets 

By measurements of the Mossbauer effect in iron, Hanna 
et al. (22) have shown that the ground state (/ = l / 2 ) of the Fe 57 
nucleus in ferromagnetic iron is split at room temperature by 
about 45 me due to the interaction between the small nuclear 
moment and the large (^330,000 gauss) magnetic field at the 
nuclear site. Following this observation Gossard, Portis, and 
Sandle (23) and Robert and Winter (24) observed this splitting 
directly by detecting the nuclear resonance of the Fe 57 nucleus 
in the effective internal field (H) . 

In a pioneering paper, W. Marshall (25) presented a theory 
for the origin of the field at the nucleus. All the sources of the 
field in his theory fluctuate very rapidly compared to the Larmor 
frequency. Thus the nucleus responds to the time-average field 
from each source. In the case of iron, which is cubic, the follow- 
ing effects contribute to the total field at the nucleus: 

(1) The demagnetizing field DM, the Lorentz field, and the external 
field HQ. M is the magnetization in each domain. 

(2) The field produced by hyperfine interaction with the 4s conduction 
electrons which are polarized by the ferromagnetically-aligned 3d-electrons. 

(3) The field produced by the hyperfine interaction with 4s-electrons 
admixed into the 3c?-baiid. 

(4) The field produced by the orbital angular momentum of the 3d- 
elcctrons. Normally, the orbital angular momentum of these electrons is 
quenched by the crystallize electric fields. However, the spin-orbit coupling 
permits the ion to have a nonvaiiishing component of the orbital angular 
momentum if the spin is oriented. This orbital contribution to H is pro- 
portional to the time-average spin orientation, i.e., to the magnetization 
M. 

(5) Finally, we have the contribution from the core polarization of the 
3s-, 2s-, and Is-electrons produced by the spin-aligned ^-electrons. 

The field at the nucleus has been found (26) to be in a direc- 
tion opposite to that of the magnetization vector. Since (5) is a 
large negative effect, it is generally believed that the core 
polarization is by far the dominant effect in iron. 

All the sources of the field considered by Marshall are pro- 
portional to the mean spin per atom, and hence to the mag- 
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netization (M) . Thus one concludes that the resonance frequency 
v is proportional to the saturation magnetization as 

v = AM (1-20) 

On the basis of equation (1-20) it was believed that precise 
measurements of the temperature-dependence of v which were 
possible by measurements of the nuclear resonance frequency 
could be used to obtain very accurately the temperature- 
dependence of M. Naturally it was assumed that A was a con- 
stant independent of the temperature. Some doubt as to the 
validity of this basic assumption was raised by the results of 
Robert and Winter (24), who showed that at atmospheric pres- 
sure v decreased more rapidly with temperature in the range 
to 300 K than M (27). The French workers noted that v could 
be expected to be a function of the temperature indirectly 
through the effect of the thermal expansion, and that different 
effects of the thermal expansion on M and v might be responsible 
for the apparently different temperature-dependences of M 
and P. They believed that if the temperature-dependences of 
these two quantities were corrected to constant volume, then 
v(T) and M(T) would be directly proportional. By making 
measurements of the pressure-dependence of v and M at a num- 
ber of different temperatures this supposition can be checked, 
since 



Thus, if the compressibility, thermal expansion, and pressure- 
dependence of v are known as a function of the temperature, 
equation (1-21) can be integrated to obtain (v(T)) v . A similar 
procedure holds for M( T) . 

G. Benedek and J. Armstrong (28) have measured the pressure- 
dependence of the resonance frequency v in ferromagnetic iron 
at -77 C, C, and 84.2 C in the range 1 to 10,000 kg/cm 2 . 
The results at C are shown in Figure 3. It should be noted 
that at each pressure the Fe 57 resonance frequency has been 



SOLIDS 



19 




2000 



4000 6000 8000 10,000 
P (kg/cm 2 ) 



Figure 3. The Fe 67 nuclear resonance frequency vs. pressure in ferromag- 
netic iron at 0C. 



measured to better than 1 kc, while the line width was ^-40 kc. 
It was found that (l/j>)(dj>/5P) r was nearly independent of the 
temperature and equal to 

(l/iO(di>/dP)oc = (~ L64 =*= - 01 ) X 10- 7 /(kg/cm 2 ) 

Before applying these results to equation (1-21) to determine the 
temperature-dependence of v at constant volume, we may briefly 
examine the source of the pressure-dependence of v. 

Since v is nearly independent of the temperature, we have 
from equation (1-18) that 



= (l/Ao)(&Ao/>P) + (l/M )(dM /dP) (1-22) 

where M and A Q are the values of M and A at T = K. We 
observe here a marked difference between the ferromagnetic and 
antiferromagnetic cases. In the case of MnF 2 , the pressure- 
dependence of the sublattice magnetization arises out of the 
pressure-dependence of the Neel temperature T N , and M is a 
constant nearly independent of the pressure, being given es- 
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sentially by 2@S per atom (/? is the Bohr magneton, S the spin 
quantum number). However, in the case of ferromagnetic 
metals each atom does not contribute an integral number of 
Bohr magnetons to the saturation magnetization. It is believed 
that this situation may arise out of an unequal distribution of 
electrons into spin-up and spin-down energy bands. Under these 
conditions we may expect that a change in the volume of the 
solid will alter the distribution of electrons in these bands and 
hence change Mo. In fact, J. S. Kouvel (29) and F. Galperin (30) 
have observed the pressure-dependence of M . As far as the de- 
pendence of M on the Curie temperature is concerned, we must 
remember that since T c ~ 1000 K, at 300 M is only weakly 
dependent upon T c . Of course, if the temperature were raised 
higher we could obtain information about the pressure-de- 
pendence of T c from the temperature-dependence of d?>/dP. 
However, in the case of iron, L. Patrick's (31) measurements of 
dT c /dP show that the. Curie temperature in iron is essentially 
independent of the pressure. 

We may now use measurements of dM /dP to determine 
(l/,4o)(dAo/c)P) from equation (1-22). J. Kouvel has found (29) 
that 

(l/Mo)(dM />P) = (~2.78 0.25) X 10- 7 /(kg/cm 2 ) 

Thus we find (l/4 )(d4 /dP) = (1.14 0.25) X 10~ 7 /(kg/ 
cm 2 ),ordlnyl/dlnF= -0.19. 

Returning now to equation (1-21), we may obtain (v(T)) v 
from the knowft temperature-dependence of (dp/dP), (1/V)- 
(dV/dT) and V(dP/dV). By carrying out an exactly similar 
procedure we may also obtain (M(T)) V . The results for (M(T)) V 
and (v(T)) v are shown in Figure 4. It is clear that even when 
corrected to constant volume, v and M are not directly propor- 
tional to one another, i.e. the coupling constant A is an explicit 
function of the temperature. In fact the experimental data indi- 
cate that the temperature-dependence of A is given by 

A = A Q (l - 0.77 X lO- 7 ! 72 ) (1-23) 



SOLIDS 



21 



The T 2 temperature-dependence of A suggests a band de- 
scription for the electrons responsible for the hyperfine interac- 
tion. In fact, Kittel has suggested (32) that the temperature- 
dependence of A may be due to the excitation of electrons higher 
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Figure 4. (v VK)/VV and (M Mo)/M in iron as a function of the 
temperature at constant volume, using Kouvel's data to remove the effect 
of thermal expansion on M. 

into the d band as the temperature is increased. The effective 
nuclear field is the resultant effect of the hyperfine interaction 
with all the spin-up and spin-down d electrons, each of which has 
a different energy, e, a different Bloch function, and hence a 
different coupling constant, a(e). This band description is very 
similar to the situation that obtains in the Knight shift. We may 
write the resonance frequence v as 
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where n (e) is the density of states at energy e for electrons with 
moments either up or down relative to the Weiss internal field. 
In this equation a^(e) is the coupling constant which essentially 
gives the field at the nucleus set up indirectly by a rf-electron with 
energy e and moment parallel or antiparallel to the Weiss field; 
/o is the Fermi function which gives the probability that a given 
energy state will be occupied ; f is the Fermi energy. We must of 
course include in the argument of/ the potential energy of the d- 
electron. This potential energy is written as &0V, i.e., it is taken 
as proportional to the saturation magnetization M = 0<7, where 
S = the Bohr magneton and 

N+ - N- 

cr = 



N+ + N- 



It is to be noted that the expression for a is very similar to that 
for v. In fact if the coupling constant a(e) were independent of 
eneigy we would immediately have v = aM/p. Thus in the 
present theory the temperature-dependence of A results from 
the energy-dependence of the coupling constant a(e). In general 
we may look upon equation (1-25) as an equation for a in terms 
of f and T. A second equation for o-, f , and T is given by the 
normalization condition 
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Equations (1-25) and (1-26) suffice in principle to determine <r 
and f as functions of the temperature. In fact these equations 
with /i + (e) ~ 1/2 are the basis of Stoner's "collective electron" 
model for the d electrons in a ferromagnet. The insertion of f 
as a function of <r and T as determined from equations (1-25) 
and (1-26) into equation (1-24) permits a determination of the 
dependence of v on a and the temperature. This procedure has 
been carried out explicitly (28) under a number of simplifying 
assumptions, viz. n + (e) = rc~(e), a + (e) = a~(e). The result gives 
= Aa = AM/P, with A containing a term in T 2 as the leading 
term in the power-series expansion of the temperature-de- 
pendence of A. As yet it has not been possible to reliably 
evaluate the coefficient of this T 2 term because of an absence of 
information on the energy-dependence of a(e) and n(e). Thus, 
at this point it is not possible to determine whether this theory 
is capable of accounting quantitatively for the observed tem- 
perature-dependence of the hyperfine coupling constant. 

However, the present theory does provide at least qualita- 
tively an explanation for the observed temperature-de- 
pendence of A in the following way. The observed constant is 
the average of contributions from d electrons distributed 
throughout the d band. It should be remembered though that d 
electrons at the bottom of the band and those at the top have 
very different radial distributions, and hence each of these pro- 
duce quite different fields at the nucleus. As the temperature is 
increased the tail of the Fermi distribution rises higher into the 
d band, thereby bringing in contributions from states with 
coupling constants different from the average. Since the number 
of states brought in is proportional to T 2 , the average value of A 
is altered by an amount proportional to T 2 . 

Kushida has recently measured (34) the pressure-dependence 
of the resonance frequency of the Ni 16 nucleus in ferromagnetic 
nickel at 273 K and 372 K. He found that the frequency in- 
creases very strongly on increasing the pressure, and that 

c = +9.20 X 10- 7 /(kg/cm 2 ) 
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When this result is incorporated with a determination of the 
pressure-dependence of M for nickel, it will be possible to de- 
termine the temperature-dependence of the hyperfirie coupling 
constant A in this ferromagnet as well as in iron. Kushida (35) 
has also studied the pressure-dependence of the nuclear resonance 
frequencies of the constituent nuclei in alloy systems like the 
Co 59 resonance in cobalt-rich alloys with Fe, Ni, Cu, Cr, Mn, 
Al, or the Fe 57 resonance in iron-rich alloys. Also, the Cu 65 and 
Cu 85 resonances have been observed in ferromagnetic Fe and Co 
alloys. 

In a ferromagnet the magnetic field in which the nucleus 
precesses is provided by the aligned electron spins within the 
solid. No external magnetic field is required. Under these con- 
ditions it would seem possible to carry out nuclear resonance 
experiments in the pressure range of 100 kilobars with the aid of 
recently developed "superpressure" presses. Such an experiment 
has been successfully carried out by J. D. Litster and G. B. 
Benedek (35a). The Fe 57 nuclear resonance frequency in 
natural-abundance iron was measured as a function of pressure 
to 65 kilobars. The upper pressure limit in this experiment was 
set by the strength of the press employed. The sample was iron 
powder whose particle size was ~10 to 20 microns. This powder 
was dispersed in high dilution AgCl powder, which served as 
the pressure-transmitting medium. The absence of appreciable 
line-broadening above 15 kilobars indicated that for our sample 
geometry the pressure was transmitted essentially hydro- 
statically by the AgCl. The nuclear resonance frequency (?) 
varies linearly with the pressure to 60 kilobars with 



c = -7.67 0.22 (kc/kbar) 

This result is in agreement with the data of Aim strong and 
Benedek discussed above. The linear variation of v with pres- 
sure suggests that the nuclear resonance frequency in iron can 
serve as a continuous pressure gauge at least up to 60,000 bars. 
The success of this experiment opens the way to nuclear 
resonance experiments in other ferromagnets, ferrimag- 
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nets, antiferromagnets, and those solids which exhibit pure 
quadrupole resonance. Such experiments should be particularly 
useful in detecting and elucidating the changes in electronic 
structure associated with high-pressure phase transitions. 

1.5 Paramagnetic Ions in Crystals 

Until this point we have looked upon the nucleus as the probe 
of the magnetic fields produced by the electrons in a solid. 
However the electrons themselves may be looked upon as 
probes since they have magnetic moments which interact with 
the fields set up by other electrons and nuclei. If individual 
magnetic ions are isolated by substituting them in very low 
concentration into a nonmagnetic host lattice, the ionic mag- 
netic moments will interact primarily with the magnetic field 
set up by their own nuclei and the electric field produced by the 
surrounding nonmagnetic ions. Our attention in this section 
will be devoted largely to the former interaction; the effects of 
the crystalline electric field will be treated more fully in Section 
2 of this chapter. The effect of the nucleus on the electron reso- 
nance frequency of the ion may be determined if we observe 
that the frequency (^1/T 2 ) with which the nucleus jumps back 
and forth among its magnetic sublevels is generally much lower 
than the frequency vh yp = ((AE)h yp /ft) corresponding to the 
hyperfine interaction energy. Thus the electron resonance line 
will be split into several hyperfine components, each corre- 
sponding to a particular orientation of the nuclear spin. By 
measuring the spacing of these hyperfine components the hyper- 
fine interaction constant between the nucleus and electron may 
be obtained directly. 

W. M. Walsh has studied the pressure-dependence of the 
hyperfine interaction for the *S-state ion Mn 2+ using the cubic 
crystals MgO and ZnS as host lattices. The impurity ion enters 
the lattice substitutionally for Mg and Zn respectively. The 
spin Hamiltonian which describes these ions is 

S$ + Sj) (1-27) 
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The first term describes the interaction between the magnetic 
moment of the 5 aligned electron spins (S = 5/2 for Fe 3 + and 
Mn 2+ ) and the applied field //. The second term represents the 
hyperfine interaction between the nuclear spin I and the elec- 
tron spin S. The final term describes the orienting action of the 
cubic crystalline potential on the spin angular momentum. 
This last term is roughly 10 times smaller than the hyperfine 
interaction term. We shall suppress it temporarily and consider 
the electron resonance spectrum that arises out of the first two 
terms in equation (1-27). In high magnetic fields (i.e., A/gftfl 
<<C 1) we may describe the nuclear-electron system in terms of 
the quantum numbers /, S, m/, and ms- The eigenvalues of 3C 
are readily calculated to the second order in the perturbation A 
to be 

/, m s ) = g$Hm s + Am T m s 

{m s (I(I + 1) - m?) - m I (S(S + 1) -ro|)} (1-28) 



equation (1-28) shows, to second order in A, how the Zeeman 
levels of the electrons are split into 27 + 1 levels by the hyper- 
fine interaction. If a microwave magnetic field is applied per- 
pendicular to the large field H, transitions can occur between 
levels satisfying the selection rules Ams = 1 and Am/ = 0. 
Since the spectrometer frequency V Q is normally kept constant 
while the magnetic field H is swept, resonant absorption takes 
place for those fields //(m/, ms) which satisfy the condition 

hv = E(m S9 m r ) - E(m s - 1, m/) 



g&H(m r , ms) + Am f 

(A*/2gfi{l(I + 1) - mf - m r (2m s - 1)} 
m.) 



In general the solution of this equation for //(m/, ms) involves a 
quadratic equation. However, if we keep in mind the fact that 
the last term is very small compared to the previous terms, we 
may to a good approximation replace H(mi, ms) in the de- 
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nominator of the last term by the first order solution of equation 
(1-29), i.e. by (hv Amj)/g^. Thus we have the explicit rela- 
tion for the field at which resonant absorption can occur 

= [hv Q - Am z 



The last term above shows that each absorption corresponding 
to a given m/ is further split into (2*8 + 1) absorption lines in 
second order. Actually, the cubic crystalline field term of mag- 
nitude a also splits each hyperfine line into (2S + 1) levels. The 
method for carrying out the detailed identification of each ab- 
sorption line in terms of (m/, ms) is given in Mattarese arid 
Kikuchi (37), who include the effect of the crystal field term. 
Once the various lines have been labelled in terms of ms and m/, 
the values of g, A, and a can be obtained immediately from the 
values of field at which absorption occurs. 

By measuring the effect of hydrostatic pressure on the 
paramagnetic resonance spectrum of Mii 2+ in MgO and ZnS, 
W. M. Walsh (36) has determined the pressure-dependence of 
the hyperfine coupling constant. In MgO he found that A 
-81.11 X lO-^cm- 1 and (l/A)(d4/dP) r = -0.035 X 10- 6 /- 
(kg/crn 2 ). Using the compressibility of MgO to determine the 
dependence of A on the volume V of the unit cell, he found that 
(d In A/d In V) T = +0.06. On the other hand, in the case of ZnS 
the coupling constant is a much stronger function of the pres- 
sure. In this case A = -63.73 X 10~ 4 cm- 1 , (!/A)(dA/dP) = 
-0.44 X 10- 6 / (kg/cm 2 ), and d In A /d In V = +0.35. 

The hyperfine coupling between the d electrons and the nu- 
cleus results from the * 'exchange polarization" of the core s 
electrons by the spins in the d orbitals. This is, of course, the 
same mechanism that we have discussed previously. The degree 
of polarization of the s electrons depends in detail on the radial 
distribution of the d wave function (38). On compression this 
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radial distribution is altered through the interaction of the d 
electrons with the neighbor ions, and hence A changes. In ZnS 
the relatively strong pressure-dependence of A, coupled with 
the low magnitude of A, indicates a strong interaction of the Mn 
orbitals with their neighbors, as is consistent with the covalent 
nature of this crystal. 

1.6 Ferrimagnets 

The effect of hydrostatic compression on the magnetization of 
ferrirnagnets has been studied by R. V. Jones and I. P. Kaminow 
(39) using electron spin resonance techniques. To measure the 
pressure-dependence of the magnetization they placed a spheri- 
cal sample in a microwave cavity in such a way that a number of 
different precessional modes of the magnetization can be excited 
by the microwave field. By orienting the external magnetic field 
so as to excite first one mode and then another one can deter- 
mine (40,41) the magnetic moment M per unit volume M = 
(N'fA/V) where /* is the moment per formula unit. The magnetic 
moment measured in this way is, of course, the total of contribu- 
tions from all the sublattices. Using this technique, they studied 
the magnetization vs. pressure in the range from 1 to 10,000 
kg/cm 2 for yttrium iron garnet, YIG, and for magnesium ferrite. 
They also obtained indirectly a value for the pressure-de- 
pendence of the erbium sublattice magnetization in erbium iron 
garnet, ErIG. 

In YIG the yttrium ion is nonmagnetic, so that the magnetiza- 
tion comes from the iron sublattices. In this garnet they found 
that 

d In /Wd In V = -(0.57 0.03) (1-30) 

In the erbium iron garnet Jones and Kaminow measured the 
pressure-dependence of the effective gyromagnetic ratio, g e n , for 
the electrons. Under the assumption that Kittel's formula (42) 

0eff = flfFe(l + (MFe/M E r)) 
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holds, knowing the ratio Mpe/MEr and the volume-dependences 
of </Fe and MF G (these may be obtained using the results on 
YIG), it is possible to determine the volume-dependence of 
MEF. Carrying out this procedure, they found that 

d In M Er /d In V = -8.4 (1-31) 

Reasoning further that the magnetization of the erbium sub- 
lattice is proportional to the magnetization of the iron sub- 
lattice, the coupling constant being the product of \Er-Fe (the 
iron-rare earth exchange constant) and the susceptibility XEr 
of the erbium lattice, one may write 

M E r = XEr[XEr-Fe]M F e (1-32) 

Assuming that XEF is a constant, using equations (1-31) and 
(1-30) one finds 

b In A Er - Fe /d In V ~ -1 (1-33) 

Thus the volume-dependence of the exchange coupling constant 
between rare earth and iron-ion lattices can be evaluated in- 
directly. 

2. Electric Fields in Solids 

2.1 Crystalline Electric Field Gradients : Nuclear Quad - 
rtipole Resonance in Ionic, Covalent, and Metallic 
Crystals. 

Nuclei with spin / greater than V 2 can possess, in addition to a 
magnetic dipole moment, an electric quadrupole moment, Q. 
This quadrupole moment can then interact with the electric 
field gradient at the nuclear site to produce a sequence of energy 
levels. The Harniltonian describing the interaction is 



- (eQ/2I(2I - 

W,) } (2-1) 
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where 



In the above expression r N = position of nucleus and V(r) = 
electrostatic potential at the point r. 

If the nucleus in question is at a site with axial symmetry 
along the z direction, and furthermore, if the lattice vibrations 
are temporarily neglected, the spacing between the states with 
quantum numbers m z and m z 1 is given by 



hv, = (ehQ/2I(2I - l))(d 2 F/dz 2 )(3m, - 3 / 2 ) (2-2) 

Thus the level spacings are proportional to the product of the 
nuclear quadrupole moment Q and the principle value <j> zz = 
d 2 V/dz 2 = </o of the field gradient tensor. Transitions between 
these levels may be induced by the interaction between the 
nuclear magnetic moment and an applied radio frequency mag- 
netic field. A measurement of pure quadrupole resonance fre- 
quency thus provides .information on the electric field gradient 
at the nuclear position in the solid. 

It is obvious that the application of hydrostatic pressure will 
alter the field gradients at the nucleus by changing the inter- 
atomic separations in the solid. Thus by 1955 a number of 
workers (43-47) had observed the effect of hydrostatic com- 
pression on the pure quadrupole resonance frequency. 

In order to analyze these experiments adequately it is neces- 
sary first to include the effect of the lattice vibrations. Following 
the discovery of pure quadrupole resonance, Bayer (47) found 
that the resonance frequency v was a strong function of the 
temperature, and he presented a simple theory to explain this. 
In the presence of lattice vibrations the components of the field 
gradient tensor fa in a fixed coordinate system will fluctuate. 
This fluctuation is at a higher frequency than that correspond- 
ing to the spacing of the energy levels of the system, so that 
only the average value of these fluctuating components fa will 
enter into a determination of the energy. Bayer considered a 
simple model of the vibration in which the field gradient tensor 
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components did not change in magnitude, but the tensor as a 
whole oscillated in an angle 6 around some equilibrium direction. 
In such a model the value of q = (5 2 F/dz 2 ) is reduced due to the 
averaging process to g (l 30 2 /2) (assuming <3C ?r/2). In 
Bayer's picture the observed decrease in v as T increases was 
due to the increase in the amplitude of the thermal vibrations. 
The Bayer theory is successful in qualitatively explaining the 
observed temperature variation of v but fails quantitatively. 
Kushida later generalized (48) Bayer's theory to include the 
effects of all the modes of vibration of the lattice and was able 
to show that the resonance frequency could be written as 



v = o(i - ! ZttT/e,)) (2-3) 

* 

where 4? is the amplitude of the ith normal mode of lattice 
vibration arid (1/60 is essentially a coefficient which describes 
the effect that a particular mode of vibration has on the mag- 
nitude and orientation of the field gradient tensor. For simple 
geometries 6 t can be interpreted as the moment of inertia of the 
particular mode of motion. By recognizing that the amplitude of 
the ith normal mode is determined by the mean energy of the 
corresponding oscillator, i.e. 

V** 2 ^ = WV2 + V(exp(WfeT) ~ 1)1 (2-4) 

it is again possible to predict the temperature variation of v. 
But even Kushida' s generalization of the Bayer theory was in- 
capable of predicting the observed temperature variation quan- 
titatively. The reason for this is that neither theory included 
the effect of thermal expansion on the magnitude of the com- 
ponents of y . 

To analyze this situation somewhat more carefully, it is useful 
in the sum in equation (2-3) to split the lattice vibrations into 
low- arid high-frequency terms depending whether they are 
above or below co^ = feTV/h, where T M is the lowest tempera- 
ture used in the experiments. The temperature variation of the 
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low frequency terms produces a temperature-dependence of the 
form (49) 

= , (i + bT + c/T) (2-5a) 

b = (-3&/2) E (l/e<w?) = (-3/2)M<l/co 2 e) (b) 



c = (-hV8*)Z(l/6i) = (-hV8*)M<l/9> (c) 

1 

Due to the appearance of the l/w< factor in the terms in equa- 
tion (2-4), the high-frequency terms generally contribute much 
less than the low-frequency terms to the temperature-de- 
pendence of *>, and can often be neglected in practice. From the 
definition of 6 and c it is seen that they contain information on 
the vibration frequencies and effective moments of inertia of 
the various vibration modes of the molecules in the crystal. 
VQ contains information as to the magnitude of the principal axis 
of the field gradient tensor 5 2 V/dz 2 . The failure of previous 
attempts to fit equation (2-5a) to the experimental data was 
that PO, 6, and c were considered independent of volume, that is, 
that they were constants independent of the volume of the 
solid. In fact, of course, they are not. One can estimate roughly 
using (5? /"o) ~ (8V/ V) that between < T < 300 K the fre- 
quency should change ^2% due to the volume effect alone. 
Also if *>o, 6, and c were independent of volume, the pressure- 
dependence of the resonance frequency would be zero. Actually, 
though, looking at Figure 5, which shows measurements of v vs. 
pressure for Cu20 (49), we see that this is far from being the 
case. To apply the theory of Kushida correctly to the experi- 
mental data it is necessary (49) to plot v vs. volume isotherms 
for a range of temperatures. It is presumed that V as a function 
of P for various temperatures is known. Then at each volume 
VQ, 6, and c can be determined. To check whether the tempera- 
ture variation at P = can be explained one must determine the 
volume at each temperature for which P = and put into equa- 
tion (2-5a) the corresponding values of *>o, 6, and c. If this proce- 
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Figure 5. The pure quadrupole resonance frequency of the Cu 63 nucleus in 
Cii20 as a function of pressure for three temperatures. 



dure is carried out it is possible to explain the temperature varia- 
tion of the resonance frequency quantitatively. 

But notice that here we have also obtained the volume-de- 
pendence of the internal field gradients from v Q (V). Also from 
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b(V) and c(V) we have the volume-dependence of the vibration 
frequencies and their magnitudes. This is rather detailed in- 
formation which casts considerable light on the internal electric 
fields and the lattice vibrations. 

The analysis described above has been carried out (49) for 3 
solids, Cu 2 0, KC10 3 , and para-dichlorobenzene, in experiments 
to 10,000 atm using temperatures as low as 77 C and as high 
as 100 C. The Cu 63 nucleus was observed in Cu 2 0, the Cl 35 
nucleus was studied as KC10 3 and/>-C 6 H 4 Cl 2 . 

In Cu 2 it was found that j> varies with volume roughly as 
(I/ V) (where ?( tO = ^(F/F )-- 96+ - 04 ). This is precisely what one 
would expect if the neighbors set up a Coulomb (I//* 2 ) electric 
field at the nucleus, for in such a case 



5 /1\ 1 

( ] oc 

dr VV V 



This fact confirms the belief that Cu 2 is an ionic crystal. 
Furthermore, from the volume-dependence of b and c it was 
found that 7 = < d In co f /d In V > is negative that is, 
the vibration frequencies decrease with decreasing volume. This 
result is not too surprising when it is remembered that the 
thermal expansion coefficient for Cu 2 O is negative in certain 
temperature ranges, thus suggesting a negative Gruneisen's 
constant for certain normal modes of vibration. In KC10 3 g 
changes from a (F/F )-- 022 -dependence at (V/Vo) = .98 to 
(F/Vo)-- 028 at 1.06 = V/VQ. In this substance the field gradient 
arises primarily from the C10 3 ~ covalent bond, with a few per- 
cent contribution from the potassium ions. The weakness of the 
observed volume-dependence of go is consistent with the view 
that the compression of the crystal changes the K + and C10 3 ~ 
distances without appreciably affecting the dimensions of the 
C10 3 ~ group. In KC10 3 y has the reasonable value of about 1.8 
(1.4 at V/Vo = .98 and 2.2 at V/V Q = 1.06). The experimental 
results on para-dichlorobenzene are so unusual that they are 
given in Figures 6, 7, and 8. In the same specimen these data 
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Figure 6. The pure quadrupole resonance frequency of the Cl 36 nucleus 
in para-dichlorobenzene as a function of pressure at T = 24.8 C. The a, , 
and 7 phases are polymorphic modifications which are produced by hydro- 
static pressure. 

contain two phenomena dear to the heart of the high-pressure 
physicist: the presence of polymorphic transitions, and a maxi- 
mum in the v vs. pressure curve. At room temperature and 
atmospheric pressure the a phase is observed, its frequency in- 
creasing with pressure until at 1600 kg/cm 2 the line quickly 
disappeared. The line could be found again 436 kc/sec higher 
with an intensity and line width about the same as the a phase. 
This new phase can be identified as a 7 phase (50). On increasing 
the pressure the v vs. pressure curve flattens until at 8500 kg/cm 2 
there is a suggestion that dv/dP may become negative. On de- 
creasing the pressure the 7 phase can be followed below 1600 
kg/cm 2 until at about 750 kg/cm 2 it disappears, to reappear 
once again in a /3 phase (50) 27 kc/sec below the a phase. The /3 
phase is metastable relative to the a phase at room tempera- 
ture, returning to it in a time which is very sensitive to the 
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Figure 7. The pure quadrupole resonance frequency of the Cl 36 nucleus in 
para-diehlorobenzene as a function of pressure at T = 0C. 

pressure (49). At lower temperatures we observe very clearly 
that the 7-phase v vs. pressure curve does possess a maximum. 
The maximum which we have observed in the 7 phase can be 
understood as follows: With a decrease in volume, the ampli- 
tude of vibration decreases, that is 6* decreases, thus tending to 
produce an increase on the resonance frequency. On the other 
hand, at low pressures VQ is roughly independent of volume. 
This is because the field gradient at the Cl 36 nucleus is produced 
by the C-C1 covalent bonding. Since in the Oth approximation, 
the application of pressure to a molecular crystal changes only 
the intermolecular distances, go can be expected to be inde- 
pendent of volume at the lowest pressures. However, as the 
volume gets smaller, the covalent bonding decreases by inter- 
molecular hybridization with neighboring molecules. This re- 
duces the field gradient, thereby tending to reduce the resonance 
frequency. The competition between these two effects causes the 
maximum in v. The maximum shifts to lower pressures as the 
temperature is dropped because 2 starts small at low tempera- 
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Figure 8. The pure quadrupole resonance frequency of the Cl 35 nucleus 
in the y phase of para-dichlorohenzene as a function of pressure at 77C. 

tures and it is therefore easier to overcome. Quantitatively, it is 
observed that at (V/V ) = 1.04 q varies as F, while at (V/V ) 
= 0.94 q varies as y- 036 . y is about 1 and is roughly independ- 
ent of volume. 

Recently T. Fuke, working in Professor Koi's laboratory at 
Tokushima University, has observed (51) the pure quadrupole 
resonance spectra of the I 127 nucleus in SnI 4 and the Br 81 
nucleus in /)-C 6 H 4 Br 2 under pressure up to 10,000 kg/cm 2 . In 
both these crystals Fuke found results quite similar to those 
mentioned above for /}-C 6 H 4 Cl2. The bromine and iodine reso- 
nances occur at about 200 me and 116 me respectively. Never- 
theless he found it possible to simply adapt the conical electrical 
lead for operation in this frequency region. Fuke also found that 
petroleum either could be used as a pressure-transmitting 
medium to about 1,000 kg/crn 2 at liquid-oxygen temperatures! 
This result, combined with the previous observation (49) that a 
50-50 mixture of n-pentane and 2-methylbutane will not freeze 
below 10,000 kg/cm 2 at 196 K, means that a high-pressure 
press designed to use a liquid pressure-transmitting medium 
can be used effectively as low in temperature as 90 K. Fuke 
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also studied the As 76 resonance in AssOe and the Br 81 resonance 
in KBr0 3 , and analyzed all his results along the lines described 
above. 

The negative slopes (dv/dP) which we have noted for the y 
phase in p-Cttltdz are about ten times smaller than those ob- 
served for the Cl 35 resonance in HgCL (52) and the Br 81 reso- 
nance in TiBr 4 (53). In this latter compound, Barnes and 
Engardt (53) have concluded that this negative pressure- 
dependence of the quadrupole resonance frequency can explain 
the strange temperature-dependence of the resonance frequency. 
In TiBr 4 the resonance frequency increases with increasing tem- 
perature from 250 C to 50 C, until at this temperature it 
starts to behave normally and decreases with increasing tem- 
perature. They argue that the thermal expansion causes the 
resonance frequency to increase as the temperature is raised, 
because (dv/dV) = (dv/dP)(dP/dV) > 0. This product is 
positive because both dv/dP and dP/dV are negative. However, 
the lattice vibrations tend to lower the resonance frequency. 
This effect dominates above 50 C, while below this tempera- 
ture the thermal expansion effect dominates. 

By measuring the pressure-dependence of the quadrupole 
splitting of the nuclear magnetic resonance lines of Na 23 in 
NaC10 3 and NaBr0 3 , Gutowsky and Williams (54) were able to 
show that the electric field gradient at the sodium nuclei varies 
approximately as V~ 2 . It was pointed out in this work, and in a 
later study of NaN0 3 (55), that in order to analyze this result it 
was important to keep in mind that the distortions in the unit 
cell will not in general be specified by a single parameter, V t the 
volume of the unit cell, even though the compression is hydro- 
static. They realized that in order to calculate theoretically the 
pressure-dependence of field gradient it is necessary to know 
how the pressure affects the internal parameters which describe 
the position of each of the atoms within the unit cell. This point 
was also made by W. M. Walsh (56) in his study of the effect of 
pressure on the crystalline field splitting of the ground state of 
Ni 2+ in NiScFe'HW). These observations made it clear that con- 
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cornitant measurements of the pressure-dependence of the 
x-ray structure factors may make a valuable contribution to 
studies of the pressure-dependence of the pure quadrupole 
spectra. 

The effect of pressure on the quadrupole resonance in a metal 
has been studied in the case of gallium (57). Despite the com- 
plexity of the crystal structure of this metal, unpublished 
analysis of the data indicates that the field gradient at the 
gallium nucleus varies roughly as I/V. 

2.2 Higher Moments of the Crystalline Fields 

When a free ion is substituted into a crystalline lattice, its 
electrons experience the electrostatic field set up by the neigh- 
boring ions in the lattice. The angular and radial dependence of 
the electrostatic potential which the central ion experiences 
can be described in terms of an expansion in powers of the 
quantity r//?, where r = l x x + l y y + \ z z is the distance from 
the central ion to the point at which the potential is evaluated, 
and R is the distance between the central ion and the nearest- 
neighbor ions. In the case that the neighbor ions are arranged on 
the 6 corners of an octahedron around the central ion, the lowest 
nonvanishing terms in the expansion of the potential apart from 
the spherically symmetrical term is 

Foctahedrai = D(x* + y 4 + r 4 - 3r 4 /5) + 0(r/W) + . . (2-6) 

D is a coefficient which depends on H and the detailed charge 
distribution of the neighbors. If the neighbors are regarded as 
points with charge ze then 

D = (35/4) (ze/fi 5 ) (2-7) 

It will be remembered that a similar expansion of the electric 
field occurs in the case of the interaction between the nuclear 
quadrupole moment and the crystalline electric field. However 
in that case the most important terms in the expansion are those 
corresponding to derivatives of the form (d 2 y/dr 2 ) r=5B0 . In the 
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case of point charges this leads to quadrupole coupling constants 
which vary with R as I//? 3 . In the paramagnetic resonance case 
the important quantity is the dependence of D, or (dF 4 /^ r4 )r=o, 
on R. The simple point charge approximation equation (2-7) 
leads then to a I//? 5 dependence for this quantity. 

As a result of this crystalline potential (equation (2-6)), the 
spatial degeneracy of the orbital states of the free ion is lifted, 
i.e. the atomic term values are split. If the degeneracy of the 
ground state is completely removed by the crystalline field, the 
orbital angular momentum will be quenched, i.e. the expectation 
value of any component of L will be zero. Thus as far as the 
magnetic properties of the ion are concerned the orbital con- 
tribution to the magnetic moment in first order will be zero. 
Physically what occurs is that the crystalline field exerts torques 
on the orbit of the electrons which continuously reorient the 
plane of rotation. As a result, on the average the orbital con- 
tribution to the magnetic moment is zero. Thus in first approxi- 
mation one would expect 'that the interaction between the ion 
and an external magnetic field would be described solely by the 
spin Hamiltonian 



(2-8) 

where g<> = 2.0023. A paramagnetic resonance investigation of 
the energy levels of this Hamiltonian would yield only g Q . 
However, in the next highest approximation it is easy to see that 
the application of a magnetic field slightly "unquenches" the 
orbital angular momentum through the spin orbit coupling. 
Physically, the application of the magnetic field orients the spin 
moment; the spin moment in turn is coupled to the orbital 
moment via the spin orbit interaction XL-S. Thus the oriented 
spin will partly orient the orbital moment. According to second- 
order perturbation theory (58) the magnitude of the orbital con- 
tribution ft (L) can be expected to be the order of 0(X/A), where 
X is the spin orbit coupling constant and A is the crystalline field 
splitting of the orbital levels. As a result of this effect the mag- 
nitude of the g value in the spin Hamiltonian will depart from 
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2.0023, and in general will depend upon the direction of the 
magnetic field relative to the crystalline axes. The deviation of g 
from 0o provides information jointly about the spin orbit 
coupling constant and the crystalline field splitting of the free- 
ion states. 

W. M. Walsh (36) has taken advantage of this fact to use the 
paramagnetic resonance method to study the effect of pressure 
on the crystalline field splitting. He studied the ions Ni 2 + and 
Cr 3 + substituted into the cubic crystal MgO. In this crystal the 
impurity ion is surrounded by a nearest-neighbor octahedron of 
2 ~ ions. Both Ni 2+ and Cr 3 + are in orbital F states in the free 
ion, the spin being S = 3/2 for Cr 3+ and S = 1 for Ni 2 +. 
The octahedral crystalline field splits the orbital F states into 
two triplets and a singlet, the singlet lying lowest. Denoting the 
singlet-triplet splitting by A, and using the point charge ap- 
proximation for the potential (equations (2-6) and (2-7)), it 
follows that A is given by 

A - (10/3) (eKr 4 )//? 5 ) (2-9) 

where (r 4 ) is the average value of r 4 for electrons in 3rf-orbitals on 
the metal ion. Because of the cubic symmetry of the crystal 
field, the g value is isotropic and the spin Hamiltonian is 

5C = gflS-H (2-10a) 

where 



(b) 

By observing the paramagnetic resonance spectrum under hy- 
drostatic pressure, Walsh found (36) that for Cr 2+ in MgO g 
changes from 1.97970 =t 0.00003 to 1.97995 =fc 0.00003 as the 
pressure changes from 1 to 10,000 kg/cm 2 (see Figure 9). This 
result demonstrates clearly the precision required in these ex- 
periments. The high accuracy was achieved by means of a feed- 
back loop from the output of the paramagnetic resonance spec- 
trometer to the current control on the magnet. Using this loop 
the magnetic field was locked to the center of the resonance line. 
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Figure 9. The pressure-dependence of the (/-value for the paramagnetic 
resonance spectrum of Cr 3+ in MgO. 



Any shifts in the magnetic field corresponding to the center of 
the line could be determined using a nuclear resonance probe to 
a few parts in 10 5 . Walsh found that 

V(<7 - 



= 1.14 X 10~ 6 /(kg/cm 2 ) 
1.08 X 10- 6 /(kg/cm 2 ) 



Cr 3 + in MgO 
Ni 2+ in MgO 



(2-11) 



In order to transform this into a form which shows the de- 
pendence of (g g ) on the distance between the metal ion and 
the oxygen neighbors one must know the compressibility in the 
neighborhood of the metal ion. In the case of Ni 2 + it is reasonable 
to use the compressibility of bulk MgO, since NiO and MgO 
have nearly identical lattice constants. While this is not the 
case for Cr 3+ ion, the compressibility of MgO was used there 
too to determine the dependence of (X/A) on B. Walsh found 



- g Q )/dlnR = 



(6.0 0.6) for Cr 3+ in MgO 
(5.7 0.9) for Ni 2 + in MgO 



(2-12) 
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If we look back at equations (2-9) and (2-10), and if X and (r 4 ) 
are regarded as independent of 73, we find at once that 

d In (g - 0o)/b In fl = -d In A/b In R = +5.0 (2-13) 

Thus, within the restrictions imposed by the assumptions made 
above, the simple point-charge picture seems nearly capable of 
accounting for the observed dependence of the orbital moment 
on the interatomic distances. However, an examination (36) of 
the magnitude of X/A indicates that the point-charge picture is 
incapable of yielding the observed magnitude of (g g Q ). 

In case the free ion has zero orbital angular momentum in its 
ground state (e.g., Mn 2+ and Fe 3+ ) one would think that there 
would be no interaction between the spin and the crystalline 
electric field. However, such an interaction, while small, can be 
observed in the paramagnetic resonance spectra of these S-state 
ions (59) in MgO. If one also includes the effect of the hyperfine 
interaction (AI*S) between nuclei and electrons, the spin 
Hamiltonian describing the spectrum is given by (36) 

X = 08H-S + (l/6)a(S/ + S/ + S 2 4 ) + AI- S (2-14) 

in which g is very close to the free-electron value. The term of 
present interest is the second on the right hand side of equation 
(2-14) whose magnitude is given by a, the so-called ground-state 
splitting parameter. \ detailed mechanism for this coupling be- 
tween the spin and the crystal field is not yet clearly established 
(60,61), and a physical description of the proposed mechanism is 
difficult as the calculation involves simultaneous inclusion of the 
spin orbit, magnetic spin-spin and crystal field interactions, 
carried out to the sixth order of perturbation theory. Different 
predictions (60,61) have been made for the dependence of a on 
the interatomic distances. Watanabe finds (60) aa(Dq) 2 , i.e. 
(d In a/d hi R) = 10 if the point-charge model is used. On the 
other hand, Powell, Gabriel, and Johnston (61.) have calculated 
the dependence of a on the crystalline field splitting parameter 
Dq where q = 4e (r 4 )/105 and D is given in the point-charge 
model by equation (2-7). They find (5 In a/d In Dq) = 2.1 and 
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3.5 respectively for Mn 2 + in ZnS and MgO. If we use Walsh's 
results described above for Cr 3+ and Ni 2+ in MgO to estimate 
the dependence of Dq on R for both ZnS and MgO we find 
d In Dq/d In R = 6. (This procedure is clearly questionable 
in the ZnS case.) One is led, using Powell, Gabriel, and John- 
ston's calculation for d In a/d In Dq to 



eor (2-15) 

21 Mn 2+ m MgO 

Walsh has also measured (36) the paramagnetic resonance 
spectra of Mn 2 + in ZnS and MgO as a function of pressure and 
has determined the pressure-dependence of the ground-state 
splitting parameter. Assuming that the local compressibility is 
identical to that of the host lattice, Walsh found 

^i^/^i m ~ 9 - 6 Mn 2 +inZnS /01 . N 

(5 In a/5 In *) expt = (2-16) 



which supports the Powell-Gabriel-Johnston calculation. As 
this discussion shows, these high-pressure measurements are 
very useful in establishing the correct theoretical description of 
such a subtle quantity as the ground-state splitting parameter of 
these S-state ions. 

2.3 The Isomer Shift 

The work of R. Mossbauer (62) has shown that an appreciable 
fraction of the 7-rays emitted by nuclei bound in solids are 
emitted without recoil of the parent nucleus. The solid as a 
whole may take up the recoil momentum, provided that the 
emitted 7-ray energy is not too large. Since the lifetime (T) for 
emission for the 7-ray can be quite long, the energy of the 
emitted 7-ray can be extremely well defined. For example, in 
the case of the decay of Fe 57 from its excited / = 3 / 2 state to the 
/ = l / 2 ground state r ~ 1.4 X 10 ~ 7 , hence the uncertainty in 
7-ray energy h/r divided by the 7-ray energy E = 14.4 kev is 
&E/E = 4 X 10~ 13 . The difference in 7-ray energy between 
nuclei in a source and an absorber can be observed by detecting 
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the resonant absorption of 7-rays by an absorbing foil. The 
resonance curve can be swept out by changing the velocity of the 
source relative to that of the absorber (63). In this way the 
Doppler shift is used to sweep the energy of the 7-ray. 

Using this method it has been found that the energy of the 
7-ray depends not only on the spacing of the energy levels of the 
bare, static nucleus, but also upon (a) the electrostatic interac- 
tion between the nuclear charges and the electronic charge in- 
side the nucleus (64,65,66) : the isomer shift, (6) the magnetic 
hyperfine interaction between the nucleus and its surrounding 
electrons (67), (c) the quadrupole interaction between the 
nucleus and the crystalline electric field gradients, and (d) the 
temperature. This last dependence comes about through the 
time dilatation produced in second order by the vibration of the 
nuclei about their equilibrium positions (68). 

The Mossbauer effect has proven to be a most useful tool for 
the study of solids because the 7-ray energy carries with it in- 
formation on the electronic structure and atomic motion in the 
solid. Furthermore, the sharpness of the 7-ray emission line 
enables a measurement of exceptionally small changes in the 
energy of the 7-ray. In fact, R. Pound, G. Benedek, and R. 
Drever (69) have measured the effect of hydrostatic pressure on 
the energy of the 14.4-kev 7-ray from ferromagnetic iron. The 
7-rays emerge from the high-pressure vessel through a conical 
beryllium window. While this experiment is not a magnetic 
resonance experiment, it provides information so like that which 
we have discussed previously that it seems appropriate to in- 
clude it. 

Since ferromagnetic iron is cubic, the quadrupole interaction 
in the / = 3 /2 state is absent. Furthermore, the method of 
observing the spectrum is to superpose the hyperfine pattern in 
the source over that in the absorber so that in fact one detects 
changes in the center of gravity of the hyperfine pattern of the 
source. Thus, changes in the energy of the 7-ray come solely 
from changes in the quantity : 

kv = &?o + AviBom + kv r el (2-17) 
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The term hv^ represents the spacing of the energy levels of a 
stationary nucleus free of encircling electrons and external 
electric and magnetic fields. The /ij>i 80m represents the electro- 
static interaction between the nuclear charge and the electronic 
charge density inside the nucleus. This interaction depends on 
the effective nuclear radius, so that if this radius is different in 
the excited and ground states there will be a change in the elec- 
trostatic interaction energy which appears in the energy of the 
7-ray. To get some feeling for this effect one may consider a 
spherical uniformly charged nucleus (charge Ze) interacting 
with the electronic charge density at the nucleus ^(O)! 2 . The 
difference in eneigy between the excited and ground state is 



where p ex and p gr are the radii of the excited and ground states 
respectively. Of course in general we have 



where the constant K can be evaluated (66) by observing the 
difference in isorner shifts between iron atoms in different de- 
grees of ionization. If one then employs the free-ion wave func- 
tions (70) to calculate the difference in |^(0)| 2 , K can be found. 
The size of hv iso m relative to hv$ is as follows: The core electron 
contribution (66) in iron is (j>isomAo)core = 2 X 10~ 8 ; the 
outermost 45 electrons give (^i s om/^o) = 5 X 30~ 12 . The minus 
sign comes from the fact that the "radius" of the excited state in 
iron is smaller than the ground state (66). While these are very 
small contributions, it should be kept in mind that the 7-ray 
emission line is extremely narrow; our estimate above shows 
&E/E~ 4 X 10 - 13 . 

The relativistic contribution due to the time dilatation pro- 
duced by the oscillation of the Fe 57 atom about its equilibrium 
lattice position is calculated as follows: 

frel = * (1 - <U 2 )/C 2 ) 1/2 ~ ?o (2-18a) 

^ - p (M(v*)/2Mc 2 ) (b) 

^ v Q (E/2Mc*) (c) 



SOLIDS 47 

where M is the mass of the atom, and E is the average vibra- 
tional energy of an iron atom. The size of v Te \ in terms of j>o is 
(W*o) = -7X10- 13 . 

By making measurements at 1, 2000, and 3000 kg/cm 2 the 
pressure coefficient of the 7-ray energy was found to be (69) 

(\/v)v/*P) T = -(2.61 0.01) X 10~ 18 / (kg/cm 2 ) (2-19) 

Thus on applying 3000 kg/cm 2 , our maximum pressure, the 
7-ray energy changed only 7.8 parts in 10 16 ! This is by far the 
smallest effect ever measured under high pressure. The 7-ray 
energy changes under pressure through two effects. The first is 
the change produced in |^(0) | 2 by altering the volume of the unit 
cell; the second is the change in the Debye temperature, which 
in turn alters the average vibrational energy E. An estimate of 
the latter effect can be made. Using the common assumption (71) 
that (5 In 9/5 In V) = 7, where 7 is Gruneisen's constant, 
one finds at 295 K that the effect of pressure on the relativistic 
term is (69) ~ .126 X 10~ 18 /(kg/cm 2 ) which is less than 5% 
of the measured pressure coefficient given in equation (2-19). 
Thus the volume-dependence of the isomer shift must be re- 
sponsible for the observed pressure coefficient, i.e. 

2> In F)(d In F/dP) (2-19) 



Using a value of K obtained as described previously, along with 
Bridgman's value (72) for the compressibility of iron, we find, 
using atomic units (1 a.u. = 0.5292 X 10 ~ 8 cm), 

(d|^(0)| 2 /d In V) = -(2.7 0.4) a.u.~ 3 (2-20) 

We may combine the present result with that discussed pre- 
viously (Section 1.4) in the case of nuclear resonance in iron. 
In that experiment one obtains the volume-dependence of the 
difference between spin -up and spin-down charge densities. 
Using Kouvel's data for the pressure-dependence of the mag- 
netization, one finds 



= -0.19 (2-21) 
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From the magnitude and direction of the magnetic field at the 
nucleus one may deduce that 

|*t(0)| 2 ~ ki(0)| 2 = 0.63 a.u.- 3 (2-22) 

Thus 

= _ OA2aiU- _ 3 (2 . 23) 



i- 

din V 
while equation (2-20) gives 

^(0)l' + 1^0)1'.. (2 . 7 0.4) a.0.-. (2-24) 
5 In V 

Note that f denotes a spin in the same direction as the mag- 
netization of the solid. From equations (2-24) and (2-23) we 
may conclude that on changing the volume of the unit cell, the 
major effect consists in changing both the spin-up and spin- 
down charge densities by about the same amount, namely 
d|ft(0) | */& In V ~ d|irt(0)| 2 /dlnV -- 1.3 a.u.- 3 . However, the 
spin-up charge density at the nucleus changes slightly more 
strongly than the spin-down density. Quite apart from the 
individual volume-dependences of the spin-up and spin-down 
charge densities, we may estimate the magnitude of the volume- 
dependence of the total charge density at the nucleus in the 
following way. If we assume that only the 4s wave function 
|^4s(0)| 2 is changed by the volume compression, and that |^4 S (r)| 2 
scales with the volume, i.e., (d In |^ 4s (0)| 2 /d In V) = 1, then 
the estimate (66) that |^ 4s (0)| 2 ~ 3 a.u.~ 3 gives (d|^(0)| 2 /d In 
V) estimate = 3 a.u.~~ 3 . This estimate compares favorably with 
the experimental result in equation (2-20). These results show 
that the combined techniques of nuclear magnetic resonance, the 
Mossbauer effect, and high-pressure physics can provide very 
detailed information concerning the electronic structure in iron. 

3. Relaxation, Diffusion and Hindered Rotation in Solids 
3.1 Nuclear Relaxation in the Alkali Metals 

We have considered previously the information which can be 
obtained from precise measurements of the frequency positions 
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of absorption lines in the magnetic resonance spectra. Consider- 
able information, however, is contained in measurements of the 
line shapes and relaxation times. These quantities can provide 
rather detailed information on atomic and molecular motion in 
solids. 

The work of D. Holcomb and R. E. Norberg (73) is a fine ex- 
ample of the use of nuclear magnetic resonance at atmospheric 
pressure as a tool for the study of atomic diffusion in solids. 
These workers measured the temperature-dependence of the 
nuclear relaxation times of the alkali metals (lithium, sodium, 
arid rubidium) over a wide range of temperature. By careful 
analysis of their data they determined the magnitude and tem- 
perature-dependence of the self diffusion coefficient D for lithium, 
sodium, and rubidium. 

The self diffusion coefficient, D, may be defined in terms of the 
partial differential equation which governs diffusion. If P(r, f) is 
the probability that a diffusing particle may be found at time t 
at position r, then P is governed by the equation 

dP/d/ = DV 2 P 

The physical significance of D may be found by considering the 
solution of the diffusion equation subject to the initial condition 
that at t = the particle was definitely at r = r . Under these 
conditions P is given by 

P(r, = {l/(47rD/) 3 > 2 } exp {-(r - r ) 2 /4Di} 

From this we see that the mean square displacement of the par- 
ticle increases linearly with the time as (r r ) 2 = 2Dt. The 
diffusion coefficient D determines the rate at which the mean 
square displacement increases. One can provide alternatively a 
microscopic description of the diffusion coefficient by turning 
to the microscopic random walk process which leads to the dif- 
fusion equations. This analysis shows that D is related to the 
mean square step length A 2 and the correlation time T required 
for a diffusion jump by the equation 

D = cA 2 /r 
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where c is a numerical constant of the order of 0.1 to 
1.0. 

Stimulated by the work of Holcomb and Norberg, R. G. 
Barnes, R. D. Engardt, and R. A. Hultsch (74) measured the 
pressure-dependence of the nuclear resonance line width in lithium 
to 3000 atmospheres, and from this estimated the pressure-de- 
pendence of the diffusion coefficient. This experiment showed that 
large and easily detectible changes in the nuclear relaxation times 
could be produced by the application of pressure to the alkali 
metals. As a result Hultsch and Barnes (75) undertook a de- 
tailed study of the relaxation times using the methods of spin 
echoes (76,77). These workers succeeded in determining the 
pressure-dependence of the diffusion coefficients in lithium and 
sodium. 

The ideas on which this more complete work is based are as 
follows. The approach to thermal equilibrium of the nuclei can 
be described in terms of two relaxation times, 7\ and TV The 
spin-spin relaxation time (T 2 ) denotes the time required for the 
nuclear spins to establish thermal equilibrium between one 
another at some spin temperature, TV The quantity T 2 also 
characterizes the time during which phase coherence can be 
maintained between nuclear spins processing in the absence of an 
rf field. The spin-lattice relaxation time, T\, denotes the time re- 
quired for the spin system to come to the temperature of the 
lattice. In simple physical terms 1\ is the time constant for the 
establishment of the full nuclear magnetization after applying a 
d.-c. magnetic field to an initially unmagnetized nuclear spin 
system. 

Both TI and T% arise from two sources: (i) the nuclear dipole- 
dipole interaction, and (2) the hyperfine interaction between 
nuclei and the conduction electrons. Since the reciprocals of the 
relaxation times are proportional to transition probabilities, and 
since the transitions produced by mechanisms (1 ) and (2) are 
independent, it follows that 

= (I/TO. + (i/T\)d (3-i) 
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where (l/Ti) e is the contribution from the conduction electrons, 
while (l/7 T i)d comes from the dipolar interaction between 
nuclei. A similar equation holds for T 2 

(l/r a ) = (1/T 2 ). + (l/T 2 ) d (3-2) 

Since the coupling with the conduction electrons produces a field 
at the nucleus which fluctuates so rapidly that the "motional 
narrowing limit" applies (78), Ti t = T 2e . Thus, we have from 
equations (3-1) and (3-2) that 

(i/r a )d - (I/TO, = ((i/r 2 ) - (i/ro) (3-3) 

The spin-spin relaxation time, (T 2 )d, is itself constituted of two 
parts (79) 

(l/T 2 ) d = (l/T' 2 ) d + (1/270, (3-4) 

where the first term on the right hand side comes fiom frequency 
components of the perturbing dipole fields near zero frequency 
and the second term represents a broadening due to the finite 
lifetime of the nuclear spin-states. Placing equation (3-4) into 
equation (3-3) we find that 

(l/T' 2 ) d = ((1/T 2 ) - (I/TO) + (l/27\) d (3-5) 

While TI and T 2 can be measured experimentally, (T 2 )d and 
(Ti)d are the quantities which carry information on the diffusion 
constant. The theory of nuclear relaxation developed by Bloem- 
bergen, Purcell, and Pound (79) and somewhat extended to 
apply to diffusion in solids (73) shows that if the correlation time 
T C corresponding to a diffusion jump is small compared to the 
reciprocal of the rigid lattice line width Aco (i.e. if AOJT C <$C 1), 
then 

(I/T' 2 ) d = 4 7 h 2 /(/ + l)N/5Dd (3-6) 

where N = number of nuclei per unit volume, D = the self 
diffusion coefficient, and d = distance of closest approach of two 
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Figure 10. The dipolar relaxation time ( T 2 ')<* in lithium as a function of 
pressure for various temperatures. The average values of the activation 
volume per mole for the five temperatures is 3.6 rt 0.3 cm. 3 

nuclei. On the other hand, if o> is the Larmor frequency (Ti) d is 
given by 

(1/TOd = 1447T7^h/(/ + l)M)/25a> 2 rf 5 (3-7) 

provided that (cor c ) 2 )>> 1. In the motionally-narrowed limit 
(coT f ) 2 l, (T^ is given by 

(l/Ti) d = 27T7^h 2 /(/ + l)7V/5Z)rf (3-8) 

In the pressure experiments of R. A. Hultsch and R. G. 
Barnes (75), care was taken to operate in a temperature range 
in which equations (3-7) and (3-6) apply. Under these conditions 
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(7\)d is large and (l/2T r 1 ) d makes a small contribution ( 
to (l/jT 2 )d in equation (3-5). Thus measurements of the pressure- 
dependence of TI and T 2 are sufficient to determine (T^d as a 
function of pressure. Using equation (3-6), this provides the 
pressure-dependence of the diffusion constant D. 

In Figure 10 we see the results of Hultsch and Barnes for the 
pressure-dependence of (T' 2 )d in lithium to 7000 kg/cm 2 . It will 
be noted that this quantity, and hence D, decreases by a factor 
of three for 7000 kg/cm 2 . In sodium a pressure of 3000 kg/cm 2 
changes (T^ d and D by a factor of about six. 

The diffusion coefficient, J9, may be written as (80) 

D = Aza 2 "o exp (-AG a /#T) (3-9) 

where n is a constant which depends on the crystal structure, a is 
the lattice constant, and VQ is the Debye frequency which is con- 
nected to the Debye temperature 9 by hv Q = k&. AG a is the 
Gibbs free energy of activation per mole. The pressure-de- 
pendence of D results from the effect of compression on a, VQ and 
AGa. We may, in fact, relate the pressure-dependence of these 
quantities directly to that of (T' 2 )d as follows: Using equations 
(3-9) and (3-6) it follows that 

AG a //?r = 6 In a + In ? - In (T' 2 ) d + In const. (3-10) 

The last term in equation (3-10) involves quantities that are in- 
dependent of the volume of the solid. The pressure-dependence 
of T' 2 is thus connected with that of AG, a, and VQ by 



- 2(dln 
+ (d In Vd In V)(d In F/dP) y - (d In (r' 2 )d/dP) r (3-11) 



where we have used a 3 = (const.) (F). If, furthermore we use 
the definition of the compressibility & = ( d In V/dP) and, as 
in the discussion of the isomer shift, use (d In v$/~b In V) = 7 = 
Gruneisen's constant, then following Nachtrieb (81) and defining 
the activation volume AF a as 

AF a == (dAGo/dP) (3-12) 
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we find at once that the pressure-dependence of (T' 2 )d may be 
used to provide the activation volume by 

AV a /RT = -(d In (T' t ) d /*P) T - 0(2 - 7) (3-13) 

If the measurements of (d In (T' 2 )d/dP) are inserted in equation 
(3-13), and one compares the magnitude of this quantity with 
that of the last two terms assuming reasonable values of 7 (i.e., 
7 = 1.0 for Li and 1.3 for sodium) it turns out that 0(2 7) is 
less than 5% of (d In (T' 2 )d/dP). Hultsch and Barnes find the 
activation volume to be essentially independent of temperature, 
with magnitudes 

AF a = 3.6 0.3 cm 3 for lithium 

(3-14) 
AF a = 9.6 0.5 cm 3 for sodium 

Comparing these with the molar volume V mo \ (V mo i = Avo- 
gadro's number X atomic volume = 7V ft), one finds 



Fmoi = 0.28 for lithium 

(3-15) 
= 0.41 for sodium 

Hultsch and Barnes applied their results to the "dynamical 
theory" of diffusion (82). This theory predicts that the activa- 
tion energy for diffusion (A// a ) is related to the heat of fusion 
(A/f m ) and the change in volume on melting A V m by the relation 

AV d - (AH m /A// a )AF m (3-16) 

A// m is determined from the temperature-dependence of the 
diffusion constant at atmospheric pressure. Ilultsch and Barnes 
find this relation upheld in their work. It is worth noticing that 
by using equation (3-16) the activation volume may be esti- 
mated by using only data obtained at atmospheric pressure. 
The dynamical theory also provides a law of corresponding 
states, namely, that 

\nDa(T M /T) (3-17) 
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where T M is the melting temperature. This relation also was 
found valid. Finally the results of equation (3-15) may be in- 
terpreted (75) in terms of a hard-sphere vacancy model. In this 
model of diffusion the activation volume may be regarded as 
consisting of two parts, the first being 7V times the volume of the 
vacancy into which the atom moves. The second contribution is 
No times the increase in volume of the lattice when the diffusing 
atom is midway in its jump onto the vacancy. This latter con- 
tribution is nearly zero for the body centered cubic alkali metals, 
while the former would be 7V times one atomic volume if there 
were no collapsing of the neighbors inward onto the vacant site. 
The observation that the activation volume per atom is much 
less than an atomic volume suggests in this model that there is 
considerable "relaxation" of the neighbors onto the vacant site. 

3.2 Diffusion and Hindered Rotation in Solid Hydrogen 

Solid hydrogen consists of a regular lattice of molecules of 
H 2 . Because of the Pauli exclusion principle, those hydrogen 
molecules whose nuclear spins are in a symmetrical spin state 
(/ = 1) have antisymmetrical rotational wave function (i.e., 
J = 1, 3, 5 . . .). Such molecules are known as orthohydrogen. 
On the other hand, those molecules which have an antisym- 
metrical nuclear spin state (/ = 0) have even rotational wave 
functions (J = 0, 2, 4 . . .) and are known as parahydrogen. 
Because of molecular rotation orthohydrogen is considerably 
higher in energy than parahydrogen. Nevertheless the transition 
from ortho to parahydrogen and vice versa is very improbable, 
and ordinary hydrogen may be thought of as a mixture of two 
distinct molecules. Since the J = 1 state has 3 substates, the 
statistical weights of the 2 species indicate that the composition 
in metastable equilibrium should be 75% orthohydrogen and 
25% parahydrogen. The ortho-para conversion rate in the solid 
is 1% per hour initially and decreases with the time. Since only 
the ortho molecules have a net nuclear spin, these alone con- 
tribute to the nuclear resonance absorption. 
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Nuclear resonance studies in solid hydrogen at atmospheric 
pressure show the following characteristics. On increasing the 
temperature around 10 K the line width decreases suddenly, 
over a range of 1-2 K, from about 6 gauss to about 1 gauss (83). 
As in the case of line-narrowing in the alkali metals, which was 
discussed in the previous section, this effect has been explained 
(84,85) by the rapid self diffusion of the H 2 molecules. When the 
internal fields begin to fluctuate rapidly compared to the Larrnor 
frequency, the line narrows. By measuring the temperature- 
dependence of the line width in the region of narrowing, the 
activation energy for self diffusion and the correlation times can 
be obtained (84,85). A second marked change in the nuclear 
resonance line takes place near 1.5 K (83). At this temperature, 
side peaks on the line appear which grow in intensity at the ex- 
pense of the central peak as the temperature is lowered. This 
change in line shape has been explained quantitatively by F. 
Reif and E. M. Purcell (86). In their theory, the crystalline 
potential set up by the neighbor molecules splits the degeneracy 
of the Mj = 1,0 sublevels. In each of these levels the orbital 
angular momentum is quenched. Hence, of all the interaction 
terms in the Hamiltonian of the hydrogen molecule (87), only 
the magnetic dipole-dipolc interaction between the nuclei in the 
molecule contributes to the line shape. The line shape itself is 
found by properly averaging this interaction over the nonde- 
generate rotational states. Finally, the intensity of the nuclear 
resonance line as a function of time gives information on the 
ortho-para conversion rate. 

Clearly, all three of these characteristics of the magnetic 
resonance line can be expected to be functions of the inter- 
molecular distances, and hence the pressure. Recognizing this, 
G. Smith and C. F. Squire studied (88) the effect of pressure on 
the nuclear resonance line between 1.2 K and 14 K. Their 
pressure range was very modest. In fact it extended only to 
about 250 atmospheres. Nevertheless, they observed a marked 
effect on the temperature at which diffusion-dominated line- 
narrowing occurred. They determined from their data that the 
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application of only 230 atmospheres changed the activation 
energy for diffusion from 230 cal/mole (at 1 atmosphere) to 
370 cal/mole. This change in activation energy corresponds at 
10 K to a 100-fold increase in the correlation time for a dif- 
fusional jump. In this low-pressure region no effect was found on 
the temperature at which the side peaks appeared. About the 
same time as Smith and Squire were carrying out their experi- 
ments, W. D. McCormick and W. Fairbank carried out similar 
studies up to 5000 atmospheres. A preliminary account of their 
work (89) indicates that the temperature at which the hindered 
rotation transition occurs (i.e., 1.5 K at atmospheric pressure) 
increases with increasing pressure. They also indicate that they 
searched for an effect of pressure on the ortho-para conversion 
rate. It seems clear that the area of experimentation opened by 
these workers merits further and more detailed attention. 

3.3 Nuclear Relaxation in Other Solids 

Proton nuclear magnetic resonance studies (90) of the proton 
line in 2,2-dinitropropane, CH 3 C(N0 2 ) 2 , have shown that be- 
tween 196 C and 5 C the CH 3 group is constantly re- 
orienting itself around one of the crystalline axes. Above 5 C, 
however, the hindrance presented against motion through this 
axis disappears and the CH 3 group undergoes general reorienta- 
tion or tumbling. There is also evidence that appreciable trans- 
lational diffusion takes place above 5. Associated with the 
transition is a sizeable change in the volume of the solid. In the 
temperature range 5 C to 3 C the volume per gram takes 
a jump of about 0.056 cm 3 /gm (91). J. J. Billings and A. W. 
Nolle have studied the effect of hydrostatic compression on this 
transition, using the sudden decrease in the line width as an 
indicator of the transition. They find that in the range 1 to 1000 
kg/cm 2 the transition temperature is a linear function of the 
pressure. At 1000 kg/cm 2 the transition temperature is 12 C. 
By also making measurements of the thermal expansion and 
compressibility around the transition they were able to deter- 
mine the specific volume increase at the transition as a function 
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of the pressure. Their results show quite clearly that the onset of 
general reorientation is not determined solely by the establish- 
ment of a critical minimum specific volume. They found that the 
transition could occur at high pressure with much smaller 
specific volume change than occurs at atmospheric pressure, 
provided only that the temperature is increased somewhat. To 
be more specific, they found that under 1000 kg/cm 2 a tem- 
perature of over 40 C is required to produce a specific volume 
change equal to that at atmospheric pressure. However the 
transition takes place when the temperature is increased to only 
12 C. This experiment shows the usefulness of the magnetic 
resonance-high pressure studies combined with accurate thermal 
expansion and compressibility measurements in elucidating the 
details of phase transitions in the solid. 

Another study of phase transitions under pressure using 
magnetic resonance is that of Goodkind and Fairbank, who 
studied the nuclear relaxation times in solid He 3 under pressure. 
They observed the effect on TI and T 2 of crossing the a-/ft phase 
boundary. Their experiments were conducted at 1.37 K and 
2.2 K with pressures up to 120 atmospheres (92). 
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Chapter II 

Liquids and Gases at High Pressure 



The study of liquids and gases under high pressure may be 
subdivided into two general categories. The first is a study of the 
shifts and splittings of the nuclear resonance line centers around 
the values given solely by the external magnetic field, HQ. The 
shifts of the nuclear resonance frequencies are known as chem- 
ical shifts. These arise partly out of the diamagnetic circulation 
and partly out ot the paramagnetic unquenching of the orbital 
angular momentum (1,2). The splittings of the nuclear res- 
onance lines are produced by an indirect nuclear-nuclear spin 
coupling via the electrons of the molecule. Both the shifts and 
splittings of the nuclear resonance lines have been used ex- 
tensively (3,4) in the elucidation of molecular species. Never- 
theless, few studies have been made of the effect of hydrostatic 
pressure on the nuclear magnetic resonance spectra of liquids 
and gases. In Section 2.1 of this chapter we discuss the effect of 
hydrostatic pressure on the chemical shifts of aqueous solutions 
of the octahedral cobalt complexes (5,6), and the pressure- 
dependence of the chemical shift of the Xe 129 resonance in 
gaseous xenon (7,8). 

The second category under which the material of this chapter 
falls is the study of the effect of pressure on the diffusion con- 
stant and relaxation time of liquids and gases. This subject 
has received more attention than Section 2.1, and will represent 
the major part of this chapter. 
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1. Magnetic Fields in Liquids and Gases 

1.1 Chemical Shifts in Octahedral Cobalt Complexes 
and Xenon 

Soon after the discovery of the chemical shifts (9) it was 
observed (10) that aqueous solutions of compounds containing 
cobalt complexes showed large paramagnetic chemical shifts of 
the Co 59 nuclear resonance. The complexes in question consist 
of a Co 3 + ion surrounded by six neighbors arranged in the form 
of an octahedron. Examples of such complexes are [Co(CN) 6 ] 3 ~ 
and [Co(N0 2 )6] 3 ~~. A detailed analysis of the chemical shifts 
in these complexes was first given by Freeman, Murray, and 
Richards (11) and Griffith and Orgel (12). Their explanation 
may be outlined as follows. In the free ion Co 3+ contains six 
d-electrons. According to Hunds rules this ion is in the state 
6 Z). The nucleus in the free ion therefore experiences a magnetic 
field of about 10 5 gauss from the combined effects of the spin 
and orbital magnetic moments of the rf-electrons. However, 
when the ion is placed in the complex the interaction between 
the (/-electrons on the Co 3+ ion and the neighbor (ligand) ions 
produces a splitting of the degeneracy of the 5cf-orbitals on the 
Co 3 " 1 " into a triplet (lowest) and a doublet. The splitting be- 
tween the doublet and triplet is identical with the crystal 
field splitting discussed in Section 2.2 of Chapter I. When the 
Coulomb interaction between the six 3rf-electrons is taken into 
account, the doublet is split into two states labeled 1 T lg and l T^ g , 
the former being the lowest. We denote the total splitting be- 
tween the ground state and the l T\ g state by the quantity A. 
The splitting A in these complexes is larger than the coupling 
energy responsible for the parallel alignment of the electron 
spins. As a result, Hunds rule breaks down and the six electrons 
completely fill the lower triplet, producing a state with *S 2 = 
and (L) = 0, i.e., the net spin is exactly zero, and while the 
magnitude of L is not zero the expectation value of any com- 
ponent of L is zero, i.e., the orbital angular momentum is 
quenched. The application of a magnetic field slightly un- 
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quenches the orbital angular momentum. The interaction be- 
tween L and the magnetic field H admixes in a fraction / ^ 
0H/A of the l Ti p excited state into the ground state. The 
magnetic field set up by a single electron with orbital angular 
momentum hi is 

H, = 01/r (1-1) 

where /3 is the Bohr magneton: eh/2 me. Thus, the additional 
field at the nucleus produced by the unquenching action of the 
field H can be crudely estimated by setting /!/ ~l z ~ (0/f/A), 
which yields 



(1-2) 

A rigorous calculation using Ramsey's theory (1), and taking 
the six electrons into account, gives 

Atf/tf = o> = 32/3 2 <lA%(l/A) (i_ 3 ) 

The quantity a p is the 'paramagnetic contribution to the chem- 
ical shift. There is also a diamagnetic contribution, cr d , produced 
by the Larmor circulation of the spherically symmetric part of 
the electron distribution, but this is about 10 times smaller for 
these cobalt complexes than the paramagnetic term, <r p , whose 
magnitude is about 0.01. Equation (1-3) has been confirmed 
(11,14) by studies of a vs. A, where A is altered by changing the 
nature of the ligand (e.g., replacing Co(NII 3 )6 by Co(N0 2 )6). 
The magnitude of A for each complex is measured optically, 
while the chemical shift a = <r p + a d is measured using nuclear 
resonance. Excellent agreement with equation (1-3) is found in 
the sense that the quality (1/r 3 ) may be taken as a constant for 
all the ligands studied. Also the magnitude of (1/r 3 ) as deter- 
mined from equation (1-3) is quite near that found (15) from 
paramagnetic resonance studies on the Co 2 + ion. 

The demonstration of the validity of equation (1-3) provides 
the experimentalist with a new method of obtaining information 
on the crystal field splitting, A. The accuracy with which changes 
in cr, and hence in A, can be detected is estimated as follows. 
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Since o- ~ 1% at a resonance frequency (v) of 10 me/sec, 10 5 
cps out of v are due solely to the contribution of the chemical 
shift. Since the line width of the resonances are ^50 cps for 
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Figure 11. The change in the chemical shift [<r(P , T ) - cr(P, T )] of 
the Co 59 nuclear resonance frequency in aqueous solutions of 
as a function of pressure at 21.4 C and 39.4 C. 



66 MAGNETIC RESONANCE AT HIGH PRESSURE 

these complexes it is possible to detect a change in v amounting 
to about 2 cps, provided that sufficient field and frequency 
stability of the magnet and spectrometer is provided. This 
sensitivity to changes in frequency corresponds to changes in 
or of 2 X 10 ~ 6 . Thus, by studying the chemical shift one should 
be able to observe changes in A which are as small as a few 
parts in 100,000. 

With these considerations in mind, an experimental study of 
the pressure- and temperature-dependence of the chemical 
shifts of the complexes [Co(CN) 6 ] 3 ", [Co(N0 2 ) 6 ] 3 -, and [Co- 
(NH 3 ) 6 ] 3+ was undertaken (5,6). In order to achieve the stability 
discussed above a spectrometer-locking arid magnet-control 
system was constructed which provided a stability of 2 X 
10~ 7 in both the magnetic field and the spectrometer for periods 
of 6 to 7 hours. An example of the results is shown in Figure 11. 
There we plot the pressure-dependence of a in [Co(N02)e] 3 ~ 
for 21.4 C and 39.4 C for pressures up to 10,000 kg/cm 2 . 
In terms of frequency shifts, the Co 59 .resonance frequency 
changes about 3800 cps to 10,000 kg/cm, 2 when *>o itself is about 
10 me. Similar data were accumulated for [Co(CN)6] 8 ~~ at 
20.3 C, 68.5 C, and 98.1 C, and for [Co(NH 3 ) 6 ] 3+ at 19.9 C 
and 63.5 C. In each case it was found that (do-/dP)/> i atm 
was nearly independent of the temperature. 

We may explain the pressure-dependence of a by observing 
that the application of pressure alters the Co-ligand distance B 
and thereby changes A. In fact, one may estimate the de- 
pendence of A on this Co-ligand distance in the following man- 
ner. If first we regard (l/r 3 ) 3d as independent of/?, an assumption 
which is consistent with the observation that (l/r 3 ) 3tf is in- 
dependent of the nature of the ligand, then turning to equation 
(1-3) and neglecting the small diamagnetic contribution to a 
we have 

(lA)(d<7/dP) ^ (!A,)(dVdP) = -(!/A)(dA/dP) (1-4) 

For small changes in R, we may express the dependence of A 
on jRas 

A = constant//?" (1-5) 
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so that equation (1-4) becomes 

(l/eO(do/dP) = n(l/JF?)(dJR/dP) (1-6) 

The quaatity of interest is n, and its evaluation requires a knowl- 
edge of the compressibility of the complex, (l/R)(dR/dP). 
This quantity may be estimated using the vibrational spectra 
of the complex. From these spectra one can obtain the force 
constant, k, which corresponds to the "breathing mode," i.e. 
the totally symmetric vibration of the complex. This force 
constant relates the displacement, dR, of a ligand to the force 
applied to each ligand in the usual normal mode treatment of 
molecular vibration. In the present case one may regard the 
applied pressure, P, as exerting a force of PA/6 on each of the 
ligands. Here A is a not too well-defined quantity that represents 
the effective surface area of the complex. This area may be 
estimated roughly by circumscribing around the ligands a 
sphere which just touches the atomic radii of the outermost 
ligand ions. With this approach one obtains 

(l/R)(<W/dP) = -(l/R)(A/6k) = -(!/#) (A/6muV (1-7) 

In the last equation coi is the vibration frequency of the breath- 
ing mode and m is the mass of a single ligand. One may now 
combine equation (1-7) with the measured values of <r and do/dp 
to determine n from equation (1-6). The results of such a 
procedure are shown in Table II-l. In the next to last column 
we find the quantity n, which by equation (1-5) describes the 
/^-dependence of the crystal field splitting. The last column 
lists the values of n which are given by the simple model that 
the electric field set up by the ligands can be approximated by 
that set up by point changes at the ligand sites. The value n = 6 
is the result if the crystal field were produced by point electric 
dipoles. (The NH 3 complex is uncharged and possesses only an 
electric dipole moment.) The electrostatic interaction between 
the 3of-electrons which in part determines the magnitude of A 
is not large enough to produce a significant difference between 
the distance-dependence of A and lODg, the ligand field portion 
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of A. In comparing the last two columns of Table II-l it should 
be borne in mind that the point-charge model is incapable of 
providing an adequate value for the magnitude of the crystal 
Geld splitting. The observed dependence on R however is 
consistent with this simple model at least in the case of the 
CN and NH 3 ligands. It may be worth reminding the reader 
that similar conclusions were reached in the paramagnetic 
resonance studies of iron-group ions in ionic crystals (See Section 
2.2 of Chapter I). The crystal field splitting produced by the 
N0 2 ligand, however, appears to have a much stronger de- 
pendence on R than that predicted by the point charge model. 

In concluding this topic it should be stressed that the method 
employed in determining the compressibility of the complex is 
a crude one. This limits the precision with which n can be 
established. Nevertheless, within the framework of this method 
it is interesting to observe that studies of the pressure-de- 
pendence of the chemical shifts of these complexes in liquids can 
provide much the same sort of information as do optical and 
paramagnetic resonance studies on solids under pressure. 

A second example of measurements of the effect of pressure 
on the chemical shift is the work of R. Streever and H. Y. Carr 
(16) and E. Hunt and H. Y. Carr (17). While studying the 
effect of pressure on the relaxation times of Xe 129 in xenon gas, 
these workers discovered that the external field, H Q , required for 
resonance at constant spectrometer frequency decreases as the 
pressure of the gas increases. They find that the change dH G 
in the field for resonance is proportional to the density of the 
gas and the field f/ . Their measurements taken at 25 C can be 
summarized by the equation 

dH Q /H Q = 4.3 X 10~ 7 p (1-8) 

where p is the gas density measured in amagat units. (The 
amagat unit is a dirnensionless quantity whose magnitude is the 
ratio of the density of the gas at high pressure to the density at 
C and 1 atmosphere.) This relation was obtained over a 
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pressure range from 30 to 110 atmospheres, which corresponds 
to a range of p from 38 to 310 amagats. 

Carr and his coworkers make the following proposal to explain 
these results. Xenon is a spherical atom. When placed in a 
magnetic field the field at the xenon nucleus departs from the 
external field only because of the diamagnetic shielding pro- 
duced by the Larmor circulation. As a result the field at the 
nucleus, H n , is less than the applied field H by an amount 

(#o ~ H n )/H Q = 0.558 X 10~ 2 = a (1-9) 

More graphically, this signifies that their operating field of 
8060 gauss was shielded by 45 gauss due to the circulation of the 
atomic electrons around the direction of the field. 

If one now considers the collisions between xenon atoms it is 
clear that during the time of collision, r ( , this circulation will be 
somewhat interupted; a will decrease by an amount 5<r to some 
smaller value only to return to its original value after the 
collision is over. Since r as well as the time T between collisions 
are small compared to the period of the nuclear precession, the 
nucleus will see only the time average of the fluctuating dia- 
magnetic shielding field. If a' is the shielding constant in the 
presence of collisions, and if (T C /T) 1, it follows that 

cr' - a ^ <5<7(r e /r)) (1-10) 

where the average is taken over the energy distribution in the 
gas. If one crudely regards Bo- and r e as independent of the 
kinetic energy of the atoms and uses a hard sphere model for 
the collisions between atoms, then 

(1/r) = 4a 2 p(dfe7Vm) 1/2 (Ml) 

where a is the diameter of the atom, ra is the mass of the atom, 
and p is the number of atoms/cc. Thus 

dH Q /H Q = <r' - a = 5aT4a 2 (7r&r/m) 1/2 p (1-12) 
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The result in equation 1-12 suggests that it would be most 
interesting to study the temperature-dependence ofdHo/Hvp. 

While the treatment above is admittedly crude, it seems 
clear that further experiments along the lines begun by Carr 
will throw light on the details of atomic collisions and on the 
alteration in the diamagnetic shielding during such collisions. 



2. Relaxation and Diffusion in Liquids and Gases 
2.1 Liquids: The Hydrocarbons and Water 

The first nuclear resonance experiment performed under high 
pressure was a study of the nuclear relaxation times and molec- 
ular self diffusion constants of liquids. The experiment was 
proposed by E. M. Purcell and carried out by the author as his 
Ph.D. thesis (18,19). Measurements were made of the nuclear 
spin-lattice relaxation times, 7\, of the proton in the liquids 
n,-pentane, n-hexane, toluene, ethyl iodide, methyl iodide, and 
water, at room temperature and in the pressure range 1 to 
10,000 kg/cm. 2 Also, the effect of pressure on the self diffusion 
constant, D, of methyl iodide and water was determined. Both 
TI and D were measured using the spin echo techniques dis- 
covered by E. L. Hahn (21,20). 

The physical ideas which underly this study are as follows: 
As was mentioned previously in Chapter I, Section 3.1, the 
nuclear spin-lattice relaxation time 7\ is the time required for 
the nuclear moments to establish a Boltzmann distribution 
among the nuclear Zeeman levels produced by the application 
of a magnetic field, //o. The equilibrium difference in population 
between these levels results from transitions produced by 
fluctuating fields at the nucleus. In the case of protons in dia- 
magnetic liquids the fluctuating field at a given nucleus is 
produced by the magnetic dipole moment of nuclei within the 
molecule as well as those in other surrounding molecules. 
Bloembergen, Purcell, and Pound (22) have presented the theory 
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of nuclear relaxation. Their theory in the particular case of 
the water molecule gives the following expression* for T\ 



1 / 
7\ \2 



3 yh* T yW N 

~ 



where 



7 is the nuclear gyromagnetic ratio of the proton 

6 is the interproton distance in the water molecule 

T C is the correlation time for molecular reorientation. It is a measure 

of the time during which a molecule maintains a particular spatial 

orientation 

NO is the number of molecules/cc 
a is the molecular radius of the water molecule 
D is the selfdiffusion constant for water 

It has been assumed in writing equation (2-1) that r c 1/W 
where coo/27r is the nuclear Larmor frequency. This assumption 
is valid throughout the pressure range employed. The first term 
in equation (2-1) is tfie contribution to 7\ from the magnetic 
dipole-dipole interaction between the two protons within 
the H 2 O molecule. Random reorientations of the molecule 
modulate this interaction by changing the angle between the 
nuclear spin vectors and the line joining the two nuclei. r c 
characterizes the time for molecular reorientation. The second 
term in equation (2-1) comes from the fluctuating fields pro- 
duced by the diffusional motion of nuclear magnetic dipoles in 
molecules surrounding the molecule in question. Equation (2-1) 
can be generalized to apply to more complicated molecules, 
the only alteration being changes in the coefficients of T C and 
N /D. 

In order to make contact with the macroscopic properties of 
the liquid it is necessary to relate T C to some macroscopic pa- 
rameter. Bloembergen, Purcell, and Pound did this by noting 

* The expression for T\ given in Reference 22 contains a numerical error. 
This is corrected in the treatise of A. Abragam (Reference #27, pp. 300- 
302), whose results we quote. 
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that the time required for a translation al diffusion jump is 
connected to the diffusion coefficient D by 

T D --= 2a*/<)D (2-2) 

By further assuming that the correlation time T C for reorientation 
is proportional to that for translation, they could relate T C 
to D. This proportionality assumption is 

r c = Ar D (2-3) 

where A is regarded as a numerical constant of the order of 
unity. Under this assumption equation (2-1) becomes 

l/7\ = [( 7 4 h 2 a 2 A/3//) + ( 7r /5)( 7 4 riWo/a)](l/D) (2-4) 

The final step in the transition to macroscopic quantities is to 
make use of the Stokes-Ein stein relation, which connects the 
viscosity rj to D and the temperature T by 

D = feT/Smja (2-5) 

where k is Boltzmaiufs constant. Thus TI can be written in 
terms of the viscosity as 

1/T, = 8T7 4 h 2 [(a 8 A/3#) + N (ir/S)](n/kT) (2-6) 

for the H 2 molecule. For molecules containing more than two 
protons, equation (2-6) can be generalized to 

1/7!= (aA + V )(i7/*r) (2-7) 

Thus, according to this theory TI should vary inversly as the 
viscosity rj. Measurements that changed 77 by changing the tem- 
perature did in fact confirm the validity of equation (2-7). 
However, it is by no means obvious that the transition from the 
microscopic equation (2-1) to the macroscopic (2-6) will be 
valid when r c and D are altered by changing the pressure. 
Bridgman had already measured (23) the effect of pressure on 
the viscosity of a number of liquids, and had shown that t\ 
increases by a factor of 10 or 100 on increasing the pressure to 
10,000 kg/cm 2 . With this information already available it 
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Figure 12. The pressure-dependence of the relaxation time, T 1} for the 
protons in n-pentane at 27 G. 

seemed appropriate to measure 1\ vs. pressure to examine the 
validity of equation (2-1 )'. 

In order to carry out the experiment it was first necessary to 
develop a nonmagnetic pressure chamber with an electrical plug 
capable of transmission of radio-frequency power into the high- 
pressure region. Following a suggestion by Bridgman, beryllium- 
copper was used as the material for the pressure vessel. It proved 
successful from the start and has been a basic ingredient of 
magnetic resonance experiments at high pressure ever since. 
An electrical plug constructed of the same material was also 
designed (18,19) and found quite satisfactory. 

Typical results for 7\ vs. pressure may be seen in Figure 12, 
which gives the data on n-pentane, CH 3 (CH 2 ) 2 CH 3 , from 1 to 
10,000 kg/cm 2 at room temperature. The relaxation time de- 
creases by about a factor of 9 to 10,000 atmospheres. However, 
as is shown in Figure 13, which plots (l/i?(P)/(l/iKO)) as well as 
(Ti(P)/Ti(0)) on a scmilog plot, the viscosity increases by a 
factor of 40. The data of Figure 13 show conclusively that 7\ 
is not proportional simply to 1/r/, but that the product TM 
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Figure- 13. The pressure-dependence of the normalized relaxation time 
7\(P)/7\(0) and the reciprocal of the normalized viscosity rj(0)/rj(P) for 
pcntanc at 27 C. 



increases as the pressure rises. In the case of n-pentane, Tirj 
increases by a factor of 4 or 5 in 10,000 kg/cm. 2 Similar results 
were found for all the other liquids. 

In searching for the reason for the breakdown of equations 
(2-6) and (2-7), one must consider the possibility that the 
Stokes-Einstein relation (equation (2-5)) does not hold. To 
examine this possibility, the spin echo method (21) was used to 
measure the diffusion constant directly. It was convenient to do 
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this for water, for which 77 is known as a function of pressure. 
This measurement showed ( 19) that D was indeed proportional 
to (1/r?) over the entire pressure range. However TI was not 
proportional to (1/17). Furthermore D was also measured under 
pressure for methyl iodide. Viscosity data were not available 
for this liquid, but comparisons with 7\ as a function of pressure 
showed that TI decreased more slowly than did D as the pressure 
increased. This indicates that equation (2-4), which does not 
use the Stokes-Ein stein relation, also fails. 

The difficulty therefore appears to lie not in the Stokes- 
Ein stein relation, but rather in the assumption that r c is strictly 
proportional to T D , as given in equation (2-3). The experimental 
results in fact imply that A decreases markedly on increasing 
the pressure. For example, in n-pentane it must decrease by a 
factor of 5 or more on increase to 10,000 kg /cm 2 . This implies 
that a reduction of the free volume of the liquid increases 
T D much more than r c . This result is in fact quite reasonable, 
and may be understood in the following way: We may relate 
r c and T D to activation energies for reorientation and trans- 
lation by T C a exp(&E r /kT) and T D a cxp(&E t /kT), where AE T 
and &E t represent the barrier against reorientation and trans- 
lation respectively. At atmospheric pressure, the constancy of 
A on changing the temperature implies AE T = &E t . However, on 
reducing the free volume of the liquid it is entirely reasonable to 
expect that the barrier against translation (AE t ) will increase 
more than that against rotation (AE r }. In a classical picture, 
the barrier against rotation for a spherical molecule will be 
independent of the pressure, while that against translation will 
clearly increase. In such circumstances A = exp(A/? 7 AE t )/kT 
can be expected to decrease markedly, as is consistent with the 
experimental observations. 

In summary, it may be said that this experiment helped to 
clarify a delicate point in the theory of nuclear relaxation in 
liquids; in so doing it threw light on the details of molecular 
motion in liquids. Its final, and perhaps most important, result 
was to show that not only magnetic resonance experiments, but 
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other experiments (24) requiring magnetic fields can be per- 
formed at high pressure. 

A possible source of objection to this experiment was later 
investigated by A. W. Nolle and P. P. Mahendroo (25). These 
workers carefully removed from the liquids dissolved para- 
magnetic impurities, the most important of which is oxygen, by 
pumping on them, and by adding magnesium filings. The range 
of pressure in these experiments was from 1 to 1400 atmospheres. 
Comparison of their results in this narrow pressure range 
with those discussed above showed that while 7\ decreased 
somewhat more with pressure than did the unpurified material, 
it still did not decrease nearly as much as the reciprocal of 
the viscosity. 

Following the measurement of the pressure-dependence of the 
diffusion constants of water and methyl iodide (19) by the spin 
echo method, D. McCall, D. C. Douglass, and E. Anderson 
(26), using the same techniques, made an extensive study of the 
temperature- and pressure-dependence of D in n-pentane, n- 
hexane, n-heptane, n-octane, n-nonane, n-decano, 2-methyl- 
pentane, 3-methylpentane, 2,2-dimethylbutane and 2,3-di- 
methylbutane. While their beryllium-copper bomb was capable 
of pressures up to about 7,000 atmospheres, the experimental 
range actually employed was only about 600 atmospheres. 
From their measurements of D vs. temperature at atmospheric 
pressure they determined the activation energy for selfdiffusion 
at constant volume (E v ). This quantity is given by 



K v = -JR(c) In D/d(l/T))v (2-8) 

(See also Chapter I, Section 3.1). From the pressure measure- 
ments they found the pressure derivative of the activation 
energy, i.e., the activation volume, AV a . This quantity is given 

by 



AV a = -/rr(51u D/^P) T (2-9) 

Comparison of the activation volume with the molar volume V 
shows that the ratio &V a /V is between 10% and 14% for almost 
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all the liquids studied. This ratio rises to 20% only in the case of 
2,2-dimethylbutane. These workers also made the empirical 
observation that AF a /F was equal, within experimental error, 
to the ratio of the activation energy at constant volume to the 
heat of vaporization 7 vap , i.e. 

AV/y E v /E vap (2-10) 

2.2 Gases: Hydrogen, Methane, Ethylene, and Xenon 

The theory of nuclear relaxation in hydrogen was first given 
in the classic paper of Bloembergen, Purcell, and Pound (22,27). 
As has been emphasized before (Chapter I, Section 3-3), hydro- 
gen gas may be regarded as a mixture of two distinct molecules: 
ortho and parahydrogen. The ortho molecules have I = 1, 
and may have the rotational quantum numbers J 1,3,5 . . . , 
while the para molecules have 7 = and may have J = 0,2,4 
. . . Only the ortho molecules contribute to the nuclear reso- 
nance, and since the difference in energy between the J = 1 
and J = 3 states is about 850 K these molecules are largely in 
the J = 1 state provided that the temperature is not too high. 
The statistical degeneracy of the J = 1 state is 3, while that 
of the J = state is 1, so that normal hydrogen is in fact a 
mixture of 25% para molecules and 75% ortho molecules. 

In the free hydrogen molecule, the magnetic fields seen by the 
nuclei come primarily from the magnetic dipolc-dipole inter- 
action between protons 



hw'llrl, - 3(I,-n)(Irn)] (2-11) 

and from the interaction between the nuclear moments and the 
field set up by the rotation of the molecule as a whole, i.e., the 
spin-rotational interaction 

3C,_ r = hw'dx+IO-J (2-12) 

where J is the rotational angular momentum operator, Ii and 
12 are the spin operators for the two protons, and n is a unit 
vector along the line joining them. There are also fields at a 
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given nucleus produced by surrounding molecules. This inter- 
rnolecular contribution is negligible in the case of hydrogen both 
in the gaseous and liquid states. 

The m/ramolecular field fluctuates as the molecule undergoes 
collisions which produce transition between the three rotational 
states (mj = =t I). This fluctuating field is responsible for the 
nuclear relaxation. If T C is the correlation time for reorientation, 
then as Schwinger has shown (22) the spin lattice relaxation 
time may be written as (27) 

l/T, = [(4/3)(o/) 2 + (3/5)(co')K (2-13) 

The constants a/ and a/' have been measured with great accuracy 
using molecular beams. This gives 

o// T = H f = 27 gauss 

(2-14) 

a/'/27 = H" = 34 gauss 

where 7 is the gyromagnetic ratio of the proton. Using these 
values, we find 

l/7i = (2.74 X 10 12 )r c (2-15) 

It should be borne in mind that each collision between mole- 
cules does not produce a reorientation. Only the relatively weak 
anisotropic part of the intermolecular interaction can produce 
changes in the spatial orientation of the molecules. Thus r c 
is always longer than the mean time between collisions. From 
equation (2-15) we see that measurements of T\ provide a direct 
determination of r c . It is the object of the experiments discussed 
below to examine the temperature- and pressure-dependence of 
T C in an effort to obtain information on the nonspherically 
symmetric part of the interaction between hydrogen molecules. 

The relaxation time 7\ in hydrogen gas was first measured in 
1946 by Purcell, Pound, and Bloembergen (28), who used a 
sample at a pressure of 10 atmospheres at room temperatures. 
It remained for M. Bloom (29,30) and M. Lipsicas to study more 
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extensively the pressure dependence of 7\ at various temper- 
atures. Their pressure range was 150 atmospheres and the tem- 
perature range was 39 K < T < 300 K. The dependence of 7\ 
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Figure 14. The relaxation time T l for hydrogen gas vs. the density at 
37 K. Note that the arnagat = (density at pressure P and temperature T)/ 
(density at C and 1 atmosphere). 

on the pressure, or better yet, the density p, may be obtained as 
follows. If the density of the gas is low enough so that only bin- 
ary collisions are involved, it is possible to define a cross section, 
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Figure 15. 7\/P vs. temperature for hydrogen gas. 

0>, for the collisions producing molecular reorientation. In 
terms of this cross section we may write 



= p(ff r v) 



(2-16) 



where p is the number of molecules/cc ; (0>v) is the cross section 
weighted by the velocity and averaged over the velocity dis- 
tribution in the gas. In this formulation o> depends not only on 
the relative velocity of two molecules, but also on whether the 
collision is between ortho molecules or between an ortho and a 
para molecule. Thus o> as defined in equation (2-16) depends on 
the ortho-para ratio in the gas. In Figure 14 is shown the result 
of Lipsicas and Bloom for 7\ vs. p, where p is measured in units 
of the amagat: p = (density at pressure P)/(density at T = 
C, P = 1 atmosphere). The linear region below 500 amagats 
gives way to a more rapid increase above 600 amagats. From 
the slope of T\ vs. p data in the low density region they ha\e 
found the temperature-dependence of (<r r v) for normal hydrogen. 
This result is shown in Figure 15. Lipsicas and Bloom find it 
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difficult lo explain classically the observed temperature variation 
of (ff r v). In fact, in the region 100 K < T < 300 K they reason 
that o> should depend on velocity as (1/V). On weighting this 
by v and averaging over the velocity distribution this leads to a 
temperature-dependence of Ti/p which goes as T~ l/2 . This 
is in distinct disagreement with the observations. Furthermore, 
below 100 K, the kinetic energy per molecule is so low that the 
de Broglie wavelength becomes larger than the molecular 
diameter. This, they argue, should prevent the occurrence of 
the close encounters that are responsible for molecular re- 
orientation. Under these conditions a T v should fall as the 
temperature falls below 100. The observations, on the other 
hand, show a marked rise in this quantity as the temperature 
falls. 

Lipsicas and Bloom have also made a preliminary study of 
the dependence of (Ti/p) on the ortho-para concentration. 
This preliminary wock indicates that the temperature-de- 
pendence of Ti/p may depend markedly on the ortho-para 
ratio. They also measured, by spin echo methods, the diffusion 
coefficient /) for the ortho molecules in normal H 2 at 78 K to 
densities up to 500 amagats. They found that D is directly 
proportional to (1 /p), as is consistent with the relation (31) 

D = /(i?/m)(l/p) (2-17) 

where rj is the viscosity of the gas, m is the mass of the molecule, 
and/ is a numerical constant of the order of unity. 

The relaxation time T\ has also been studied as a function of 
pressure in methane, CH 4 , and ethylene, C 2 H 4 (32). These 
measurements for CH 4 , taken at 193 K and 293 K, give 

T^p = 3 X 10- 2 sec/amagat at T = 193 K 

(2-18) 
Ti/p = 2 X 10 - 2 sec/amagat at T = 293 K 

In ethylene 

= 0.25 sec/amagat at T = 293 K (2-19) 
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An equation similar to (2-15) has been given (32) for CH 4 and 
C 2 H 4 under the assumption that the spin-rotational interaction 
is small in these molecules. This assumption has been challenged 
by J. Waugh and C. S. Johnson (33) on the basis of recent 
molecular beam results. Using only the dipole-dipole field, 
Bloom et al (32) find that 

7\r c = 1.24 X 10- 11 sec 2 for CH 4 

and (2-20) 

7V C = 3.85 X 1C- 11 sec 2 for C 2 H 4 

From these equations and their measurements of T\ they find 
the ratio T/T O where r is the time between collisions. To 
estimate r they regard the molecules as hard spheres of diameter 
. The kinetic theory of gases then gives (31) 



They thus estimate for CII 4 (2-21) 

T/TC = 0.27 at T = 193 K 
T/TC = 0.15 at T = 293 K 

while for CsH 4 

T/T C = 0.66 at T = 293 K 

Summarizing this discussion of the experimental studies of 
relaxation in polyatomic gases, we observe that a theory 
capable of explaining the temperature-dependence of T\/ p for 
the low density region has not yet been produced. Neither has 
it been possible to account quantitatively for the deviations of 
TI from proportionality to p in the region of high density. 
Nevertheless it seems clear that a successful theory will throw 
light on the small anisotropic interactions between the mole- 
cules in the gas. 

We have seen that the nuclear relaxation in polyatomic 
gases is produced primarily by inframolecular magnetic fields 
that fluctuate from one value to another as a result of the 
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collisions between molecules. The situation is quite different in 
a monatornic gas. In this case the relaxation is produced by the 
mfcnnolecular fields that exist only during the collisions be- 
tween atoms (22,27). As a result, the nuclear relaxation rate 
1/Ti is proportional to the number of collisions per unit time. 
Thus if r is the mean time between collisions, we may expect 
(l/Ti)a(l/r). Furthermore, since l/r is proportional to the 
density p, we may expect in monatornic gases that the relaxa- 
tion rate is proportional to the density (\./Ti)ap. 

Using gas pressures between 30 and 110 atmospheres, the 
nuclear relaxation of the xenon 129 nucleus (I = V 2 ) has been 
studied (7,8) in natural xenon gas. This gas consists of 29% 
Xe 129 and 21.2% Xc 131 , which are the constituents with magnetic 
moments. The remaining atoms have nuclei with even mass 
numbers and hence no nuclear moment. These are Xe 132 , 
Xe 134 , Xe 136 , Xe 130 , and Xe 128 . The nuclear signal from the 
Xe 129 is weak, and the nuclear relaxation times 7\ are ex- 
tremely long (10 2 <7Y<3 X 10 3 sec), making a measurement of 
7\ quite difficult. Nevertheless, E. R. Hunt and H. Y. Carr (8) 
have shown that when sufficient care is taken to remove para- 
magnetic impurities, T\ is related to the density by 

l/7\ - (5.7 X 10- 6 )/> (2-22) 

where p is measured in amagats. This dependence on the density 
is exactly what one would expect for relaxation through inter- 
molecular fields. However, if one regards the interrnolecular 
field as arising out of the nuclear dipole-dipole interaction, the 
theory gives a value for 7\ which is three orders of magnitude 
longer than that observed. Thus some more effective relaxation 
mechanisms must be operative during atomic collisions. Such 
a mechanism has very recently been proposed by H. C. 
Torrey (34). 
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Appendix 

Experimental Methods 



The methods and materials required to carry out magnetic 
resonance experiments in the pressure range from 1 to 65 
kilobars and the temperature region from 1 K to 400 K arc 
necessarily quite varied. Nevertheless, we may roughly divide 
these methods and materials into the following categories: 
(a) pressure-transmitting media, (6) pressure multipliers, (c) 
pressure vessels for the samples and electrical plugs, and (d) 
pressure-measurement methods. 

A liquid pressure-transmitting medium which has a very wide 
useful temperature and pressure range is a 50-50 mixture of 
n-pentane and 2-methylbutane. Bridgman used this mixture in 
his experiments to 30 kilobars at room temperature. It can also 
be used (1) to transmit pressure to at least 10 kilobars without 
freezing at 196 K. Petroleum ether is a liquid very commonly 
used at room temperature in the region of 12 kilobars. Fuke has 
shown (2) that this liquid does not freeze at 90 K provided 
that the pressure is below 1000 bars. His observation suggests 
that it may be possible to find other liquids which do not solidify 
at 77 K, provided that the pressure is kept rather low, say, 
below 5 kilobars. Normally though, below 77 K helium gas is 
used as the hydrostatic pressure-transmitting medium. The 
melting curve of helium below 40 K has been obtained by 
Mills and Grilly (3) and between 35 K and 50 K by Holland, 
Huggill, and Jones (4). In Table A-l the freezing pressure, 
PF, of helium is listed as a function of temperature, T. The 
Simon melting equation applies rather well to the data, so that 
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TABLE A-l 

7 T (K) 1 4.2 10 20 30 40 50 
P f (atra) 25 140 580 1700 3300 5100 7300 

for extrapolation above 50 K one may employ the relation (4) 
(P/,/16.45) = (7V0.992) 1 - 554 - 1 

where P F is measured in atmospheres, T in K. If it is necessary 
to conduct experiments at higher pressures than the freezing 
pressure of helium it is possible to use solid hydrogen or solid 
helium. The work of Dugdale and Hulbert (5) is very encourag- 
ing as to the hydrostatic nature of the stress transmitted by the 
highly-compressible solid helium. Similarly, Jennings and 
Swenson (6) have found that solid hydrogen may be used 
satisfactorily as a transmitter of approximately hydrostatic 
pressure up to 10 kilobars through their experiments on the 
effect of pressure on the superconducting transition tempera- 
ture of a number of metals. In order to carry out magnetic 
resonance experiments into the 100 kilobar region AgCl can 
be used as a pressure-transmitting medium. This material 
apparently is capable of transmitting the pressure hydrostatically 
to a powder sample intimately dispersed in the AgCl. The 
sample geometry used in a magnetic resonance experiment to 
65 kilobars will be described in a succeeding paragraph (see 
Figure A-4) . 

The compression of the transmitting medium is produced by 
a pressure multiplier or high-pressure press. Such a press 
essentially consists of two pistons connected to one another. 
One piston has a large area, A\, the other has a small area, A 2. 
A low-pressure hydraulic pump applies a pressure PI to the 
large piston. This in turn pushes on the smaller piston. The 
pressure P 2 at the end of the narrow piston is therefore given 
by P 2 = Pi(Ai/A z ), i.e., the pressure generated by the low- 
pressure pump is multiplied by (Ai/A 2 ). The design considera- 
tions connected with these multipliers can be found in a number 
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of books and review articles (7,8,9,10). Presses suitable for 
gaseous or liquid transmitting media are readily available from 
a number of commercial firms. The pressures developed in such 
presses may be conveyed very conveniently to the vessel con- 
taining the sample by the use of hard-drawn Type 316 stainless 
steel tubing (11), o.d. Vs", i-d. 0.025". Warschauer and Paul 
(11) have found that this tubing, which is flexible and has a low 
thermal conductivity, can withstand an internal pressure near 
20,000 atmospheres before bursting. Presses suitable for use 
with solid hydrogen or helium as pressure-transmitting media 
have been described by Swenson (6,12) and Stewart (13). 
The superpresses for operation up to the 100 kilobar region are 
commercially available in the form of "belts" (14), "tetrahedra," 
and piston and cylinder combinations. 

Magnetic resonance experiments normally require the ap- 
plication of a magnetic field. Hence, nonmagnetic pressure 
vessels to contain the sample are usually necessary. It has been 
found (15) that beryllium -copper, a precepitation-hardening 
alloy of 98% Cu and 2% Be (by weight) may be used success- 
fully for this purpose. The specific alloy generally employed is 
solution-annealed Berylco No. 25, which is easily machined. 
However beryllium powder is toxic, so that adequate ventilation 
should be provided if grinding operations are performed. 
Normal machining under oil should cause no difficulties. After 
machining, the material may be hardened by heating at 325 C 
for 3 hours and then quenching in oil. This should bring the 
hardness up to Rockwell C-39 to C-41. Beryllium-copper 
pressure vessels have proven themselves at room temperatures 
to 20,000 atmospheres (16), and at or below 4 K to 10,000 
atmospheres. An article reviewing experience in fabricating 
Be-Cu pressure vessels has been published (16). 

Of these experiences, one hard-won lesson is particularly worth 
mentioning. The normally available Be-Cu rod or bar stock, 
instead of being homogeneous in the interior, occasionally 
contains bubbles of gas. This imperfection is known as "pipe" 
or "porosity." If such material is employed in the construction 
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Figure A-l. Beryllium-copper pressure vessel and electrical plug used 
with a liquid pressure-transmitting medium for nuclear resonance experi- 
ments to 10,000 bars. 
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Figure A-2. Small pressure vessel and electrical plug used for nuclear 
resonance experiments with helium as the pressure-transmitting medium 
below 77 K. 

of a pressure vessel, the gas bubbles act as stress risers and the 
bomb will crack under pressure. Ultrasonic testing before 
fabrication is often helpful in detecting such flaws. Neverthe- 
less, initial tests of a new bomb under pressure should be made 
with adequate protective shielding. Stainless steel has also 
been used as material for pressure vessels (16). It should be 
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kept in mind that in the case of pure quadrupole resonance or 
nuclear resonance in ferromagnets, antiferromagnets, and 
ferriniagnets, an external magnetic field is not necessary. In 
such experiments steel may be used for the pressure vessel. 

In order to bring electrical leads into the high-pressure region 
it has proved possible to adapt the conical seal of Bridgman for 
use in the radio-frequency and microwave region. In the previous 
figures below we see examples of pressure vessels and electrical 
plugs used in various frequency regions. Figure A-l shows a 
typical beryllium-copper high-pressure vessel (A) and electrical 
plug (C) used in nuclear magnetic resonance experiments (1) to 
10,000 atmospheres with liquid pressure transmitters. The 
rf current passes into the high-pressure region along small- 
diameter wire (J3) which connects to the cone (//). This cone 
seals against leak by pressing against a conical insulator (/) 
of artificial mica or lava. These conical seals have been used 
successfully from d.c. through 200 me (2). A flexible connection 
makes contact from the cone to the rf coil (L) inside the sample 
holder (J). The pressure transmitting liquid passes into the 
vessel through the stainless steel tubing (E) which is snugly 
threaded and sweat soldered into the plug at G. A soft copper 
plate may be soldered over the end of this connection to prevent 
the high-pressure liquid from forcing its way into and ultimately 
past the threads (G). The entire plug is sealed in the bomb with 
Bridgman's arrangement of packing washers (/)). The center 
washer is lead while the other two are Everdure. It is worth 
noticing that this 10,000 atmosphere system is quite small. 
Its o.d. is only I 3 // and its length is about 5 1 /2* r overall. It is 
used very conveniently with high- or low-temperature baths 
inserted within the pole faces of an electromagnet. A smaller 
beryllium-copper system used for operation (17) at tempera- 
tures below 70 K and for pressures below ^2500 atmospheres 
is shown in Figure A-2. To assist in sealing the cones at low 
temperatures against leakage of the gaseous pressure-trans- 
mitting medium (usually helium) it is helpful to smear 
some Dow-Corning silicon fluid No. 200 (viscosity: 12,500 
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Figure A-3. Pressure vessel and microwave plug used to 10,000 bars. 
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Figure A-4. Radio-frequency coil and sample assembly at the center of a 
belt- type press for a magnetic resonance experiment to 65,000 bars. 

centistokes) over the end of the cones. Also, in this system the 
central washer of the three sealing the plug should be gold 
rather than lead. The thermal contraction of gold is less than 
that of Be-Cu, thus a seal made at room temperature will 
maintain itself at low temperature. 

In Figure A-3 we show the microwave high-pressure system of 
Walsh and Bloernbergen (18) which operates to 10,000 at- 
mospheres. The system has two separate beryllium-copper 
plugs. The upper plug (A) is used to introduce the pressure 
fluid, while the lower is the microwave transmission plug 
which is designed to operate around 9 kmc. The microwave 
connection starts at a standard 50-ohm transmission line which 
tapers to 0.100" o.d. At the conical Be-Cu plug (D) the line 
tapers out to 0.250'' o.d. The cone seats on an artificial mica 
washer (W) which serves as insulator and pressure seal. The 
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Source 



To high pressure press 

Figure A-5. Small pressure vessel with beryllium window used in a 7-ray 
resonant absorption experiment. 



microwave cavity is press fitted on the end of this transmission 
line. The cavity (C) is coaxial with center stub X/4 long (T.E.M. 
mode). The coupling to the transmission line is capacitive and 
is adjusted by changing the gap between the center stub and the 
cone. The sample is placed at the upper end of the cavity where 
the microwave magnetic field is maximum and the electrical 
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field is minimum. The remaining volume of the cavity is filled 
with Teflon or polystyrene. The loaded Q of this cavity is 
about 500. For operation at higher frequencies it is possible to 
couple into the microwave cavity through a sapphire window. 
This method has been developed by Lawson and Smith (19) 
for use around 24 kmc. Since sapphire is transparent, their 
design may be used for optical as well as microwave experiments 
to 10,000 atmospheres. 

In Figure A-4 we show the central portion of a belt-type press 
(14) within which one can generate pressure up to the 100 
kilobar region. This type press has been used (20) to conduct 
nuclear resonance experiments at 45 me in ferromagnetic iron 
to 65 kilobars. The sample is powdered iron which is dispersed 
in low concentration in a matrix of AgCl. The rf coil passes 
out of the lava cylinder surrounding the sample, and through 
holes in the conical sealing gaskets. 

Figure A-5 shows a small beryllium-copper bomb with a 
beryllium window. This vessel was used to study the effect of 
pressure on the energy of the y-ray emitted from the excited 
state of the Fe 57 nucleus in iron (21). The beryllium window 
will permit transmission of x-rays or 7-rays at least up to 3500 
atmospheres. A test of this system to 10,000 atmospheres 
failed due to extrusion of the beryllium out the conical hole at 
about 9500 atmospheres. This failure may have been caused by 
the fact that the beryllium was originally fabricated by ex- 
trusion through dies along the same direction as that of the 
cone axis. A different orientation of the cone relative to the 
original extension axis may correct this difficulty. 

The measurement of the pressure transmitted by liquid or 
gaseous media has been greatly facilitated by Bridgman's 
observation (7) that the resistance of manganin wire varies 
linearly with pressure at least to 13,000 atmospheres. By 
calibrating a coil of manganin wire at one fixed pressure point 
and at atmospheric pressure the variation of coil resistance 
with pressure can be established and the coil used as a pressure 
gauge. The point most frequently to calibrate the manganin 
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gauge coils is the freezing pressure of mercury at C. This 
pressure (P M ) has been evaluated very precisely recently by 
D. H. Newhall, L. H. Abbot, and R. A. Dunn (22) using a very 
accurate controlled-clearance piston gauge (23) (sometimes 
called a dead-weight piston gauge). These workers find P M = 
7566.2 3.8 (bars) or 7715.6 =fc 3.8 (kg/cm 2 ). This result is 
about 1% larger than Bridgman's 1911 determination of P M 
7640 (kg/cm 2 ). 

For measurements of pressure in the ultrahigh pressure 
region where solid pressure- transmitting media are employed, 
the following fixed transition points have been established (24) : 
bismuth HI, 25,410 95 bars; thallium I-IIT, 36,690 100 
bars; barium II-III, 58,600-59,100 kilobars (25). These transi- 
tions may be detected by the change in the electrical resistance 
of wires made of these elements. 
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